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ABSTRACT 


Passive  bistatic  radar  (PBR)  and  passive  multistatic  radar  (PMR)  use  opportunistic 
transmitters  to  detect  and  locate  targets.  In  this  thesis,  a  maritime  scenario  was  modeled 
with  merchant  vessels  serving  as  multiple  non-cooperative  opportunistic  transmitters, 
while  a  frigate  warship  equipped  with  Electronic  Warfare  (EW)  and  Direction 
Finding  (DF)  receivers  takes  on  the  role  of  the  receiver  in  a  PBR/PMR  configuration.  The 
targets  are  assumed  to  be  the  generic  Formidable-class  frigate. 

A  MATLAB  model  is  developed  to  simulate  the  operating  environment  and 
passive  detection  and  location  process.  Detection  coverage  is  investigated  to  propose 
optimal  PBR/PMR  configurations  and  geometry,  while  elliptical  and  hyperbolic  target 
location  methods  are  explored  to  quantify  the  effects  of  PBR/PMR  parameters  and 
geometry  on  target  estimated  location  uncertainty. 
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I. 


INTRODUCTION 


A.  OVERVIEW 

Since  the  concept  of  radar  engineering  was  first  demonstrated  in  1904,  progress  in 
radar  technology  has  been  driven  by  growing  requirements  for  radar  performance  and 
rapidly  changing  operating  environment.  A  long  detection  range  and  wide  coverage, 
measurement  accuracy,  greater  system  capacity,  and  an  ability  to  operate  with  the 
presence  of  interference  are  some  of  the  fundamental  radar  requirements  that  have  been 
established  over  the  past  few  decades  [1],  These  radar  performance  characteristics — 
together  with  the  need  to  detect,  separate,  classify,  locate,  and  track  sources  of  emissions 
in  multi-target  environments — triggered  the  development  of  passive  radar  detection  and 
location  techniques. 

The  vast  majority  of  today’s  deployed  radar  systems  are  monostatic,  that  is,  the 
transmitting  and  receiving  antennas  are  collocated.  Despite  the  advancements  in 
antennas,  transmitters,  receivers,  and  processing  technology,  as  well  as  passive  radar 
systems,  conventional  monostatic  radar  remains  a  double-edged  sword — whereby  it 
detects  targets,  but  radar  transmission  makes  it  vulnerable  to  detection — and  may  not  be 
the  best  option  to  address  certain  operational  scenarios.  A  promising  solution  is  to  use 
multiple  radar  transmitting  and  receiving  sites  to  exploit  spatial  advantage  for 
coordinated  target  detection.  Multi-site  radars  can  be  broadly  classified  into  bistatic  radar 
and  multistatic  radar.  Bistatic  radar  is  a  radar  system  where  the  transmitter  and  receiver 
are  located  at  different  sites  [2],  Similarly,  a  multistatic  radar  system  utilizes  multiple 
spatially  separated  transmitter  and  receiver  sites  where  the  target  information  from  all 
receivers  is  fused  [1].  Passive  bistatic  or  multistatic  radar  capitalizes  on  transmitters  of 
opportunity  to  detect  and  locate  sources  of  transmission  or  targets  without  deliberate 
emissions.  The  illuminators  are  not  limited  to  radar  signals  and  include  (but  are  not 
limited  to)  analog  TV,  FM  radio,  digital  video  broadcasting  terrestrial  (DVB-T),  digital 
audio  broadcasting  (DAB),  cellular  network,  WiFi,  and  Global  Positioning  System  (GPS) 
satellite  signals  [3], 
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In  naval  operations,  targets  employing  low  radar  cross  section  (RCS)  and  radar- 
absorbent  material  (RAM)  design  methods,  coupled  with  the  use  of  highly  sensitive 
electronic  warfare  (EW)  receivers  on  warships,  have  changed  the  nature  of  the  game. 
Military  ships  are  pressured  to  limit  transmissions  to  avoid  detection  by  highly  sensitive 
EW  receivers.  Furthermore,  the  use  of  low  probability  of  intercept  (LPI)  radar  in  a 
monostatic  configuration  results  in  weak  returns  from  low-RCS  targets,  restricting 
detection  capability  and  compromising  situational  awareness.  This  thesis  explores  the  use 
of  multiple  pairs  of  passive  bistatic  radar  (PBR)  to  detect  low-RCS  targets  using 
opportunistic  transmissions  as  a  possible  solution  and  to  validate  the  findings  in  [4], 
Figure  1  shows  the  bistatic  radar  geometry  for  a  PBR  pair.  The  direct  line-of-sight  (LOS) 
distance  between  the  transmitter  (Tx)  and  receiver  (Rx)  is  known  as  the  baseline  and  is 
denoted  by  L  .  The  distance  between  the  transmitter  and  target  is  denoted  as  RT  while 

the  distance  between  the  receiver  and  the  target  is  denoted  as  Rr  .  The  bistatic  angle  [3  is 
the  angle  subtended  between  transmitter,  target,  and  receiver. 

The  proposed  PBR  approach  offers  potential  advantage  in  the  detection  of 
stealthy,  low-RCS  targets  which  are  designed  to  minimize  monostatic  radar  echoes.  RCS 
returns  of  stealthy  ships  vary  with  bistatic  angle  f3  and  can  be  sufficiently  large  at  certain 
return  directions.  PBR  capitalizes  on  this  characteristic  using  opportunistic  transmissions 
to  detect  low-RCS  targets.  Being  a  passive  system,  PBR  allows  the  receiver  to  remain 
covert,  making  it  more  resilient  to  detection  and  electronic  attack  (EA) — in  the  form  of 
jamming  and  anti-radiation  missiles  (ARMs).  The  ability  to  leverage  available 
transmission  and  to  detect  targets  passively  also  serves  to  enhance  situational  awareness, 
thereby  advancing  one’s  position  in  the  observe-orient-decide-act  (OODA)  loop  during 
an  operation  [4],  It  is  also  advantageous  to  use  multiple  transmitters  at  separate  locations 
for  detection  as  it  adds  spatial  diversity,  which  enhances  detection  accuracy  and  aids  in 
removing  clutter,  interference,  and  potential  system  errors.  Lastly,  the  PBR  system 
proposed  requires  no  additional  equipment  as  all  necessary  hardware  is  currently 
available  on  most  naval  ships. 
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Target  (Tgt) 


Figure  1.  Bistatic  radar  geometry. 


B.  HISTORY 

The  concept  of  Radio  Detection  and  Ranging  (RADAR)  was  first  demonstrated  in 
1904  by  German  inventor  Christian  Hiilsmeyer  in  a  monostatic  setup  [5].  Following  this, 
radar  experiments  in  the  United  States,  the  United  Kingdom,  France,  Italy,  Russia,  and 
Japan  were  carried  out  predominantly  with  bistatic  radar  operating  in  the  forward¬ 
scattering  configuration  [6],  However,  the  invention  of  the  radar  duplexer  in  1936 
addressed  the  issue  of  transmitter-receiver  isolation  and  broadened  the  application  of 
single-site  radar  on  aircraft,  ships,  and  mobile  ground  units.  By  the  end  of  World  War  II, 
bistatic  radar  was  displaced  by  monostatic  radar,  with  the  former  experiencing  periodic 
resurgences  [5], 

1.  First  Resurgence 

The  1950s  saw  renewed  interest  in  bistatic  radar  with  developments  in  missile  and 

satellite  detection,  location  and  tracking,  semi-active  homing  missiles,  hitchhiking,  and 

second-generation  forward-scattering  fences  [2],  During  this  period,  advancement  in 
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radar  theory  led  to  a  better  appreciation  of  concepts  on  match  filtering,  ambiguity 
functions,  statistical  theories  on  detection,  moving  target  indicator  (MTI)  radar,  and 
synthetic  aperture  radar  (SAR)  [7-9].  The  term  bistatic  radar  originated  with  K.  M. 
Siegel  and  R.  E.  Machol  in  1952  [10]. 

2.  Second  Resurgence 

The  development  of  counter-measures  to  anti-radiation  missiles  (ARMs)  and 
emitter  location-directed  artillery  in  the  1970s  resulted  in  the  second  resurgence  in 
bistatic  radar.  With  a  dual  or  multiple  site  configuration,  the  effectiveness  of  electronic 
support  measures  (ESM)  directed  attacks  can  be  reduced  by  locating  the  transmitter  away 
from  the  receiver  or  into  a  less  vulnerable  sanctuary  [5].  The  advent  of  digital  processing 
techniques  during  this  period  increased  the  processing  capability  in  MTI  operation  modes 
and  allowed  real-time  airborne  SAR  mapping  [7], 

3.  Third  Resurgence 

Research  on  bistatic  space-time  adaptive  processing  (STAP)  to  address  moving 
clutter  and  concepts  to  improve  bistatic  SAR  images  signaled  the  start  of  the  third 
resurgence.  It  was  also  during  this  period  that  passive  bistatic  radar  surfaced  as  a  possible 
counter  stealth  technique.  The  idea  for  PBR  is  to  utilize  commercial  broadcast  signals  in 
bistatic  or  multistatic  configurations  to  detect  low-RCS  targets  [5].  Since  then,  several 
PBR  systems  have  been  developed  and  evaluated  for  air  surveillance.  Some  notable  PBR 
air  surveillance  systems  are  listed  in  Table  1. 
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Table  1.  Parameters  of  significant  passive  bistatic  radar  programs  designed 
and  tested  for  air  surveillance.*  After  [5]. 


System 

Silent 

Sentry™ 

TV-Based 
Bistatic 
Radar  (I) 

TV-Based 
Bistatic 
Radar  (II) 

FM  Radio  - 
Based 
Bistatic 
Radar 

Multistatic 

HDTV-Based 

Radar 

Developer 

IBM,  now 
Lockheed 
Martin 

Univ  College 
London 

DERA,  United 
Kingdom 

NATO 

SAIC,  U.S. 
Army 

Decade 

configuration 

1980-2000 

1980 

1990 

2000 

2000 

Transmitter 

operation 

Multistatic: 

Rx:  1 

Tx:  up  to  6 

Bistatic 

Bistatic 
Near- forward 

scatter 

Bistatic 

Multistatic: 

Rx:  4 

Tx:  1 

Baseline 

100  km 
typical 

12  km 

150  km 

50  km 

10  km 
typical 

Target 

Aircraft 

Missile 

launches 

Aircraft 

Aircraft 

Aircraft 

Aircraft 
below  -  5000 
ft 

Target  data 

Range 

Doppler 

Bearing 

Range 

Bearing 

Doppler 

Bearing 

Range 

Doppler 

Bearing 

Range 

Doppler 

Measured 

performance 

Rm=  100  km 

150  km 

2- D  tracks  on 

A/C 

3- D  tracks  on 

missile 

launches 

Rm~  25  km 
Occasional 
A/C 

detections,  but 
mostly 
negative 

Rm~  160  km 
Detections  on 
high  and 
medium 
altitude  A/C 
but  only  1/3 
tracked 

Rm  ~  175  km 
Achieved  with 
innovative 
direct  path 
excision 

Rm~  30  km 
Target  location 
via  multi- 
lateration 
Ghost  excision 
via  Doppler 
association 

Status 

Version  3  for 
sale  to  U.S. 
Government 
for  <  1  million 
dollars 

Work 

continues 

Program 

ended 

Program 

ended 

Program 
continuing 
possible  in  a 
multistatic 
mode 

Test  phase 
complete 
Awaiting 
evaluation/ 
funding 

is  the  equivalent  maximum  monostatic  range  defined  as  (RrRt) —  Rm  ,  where  Rr  is  the 
receiver-to-target  range  and  RT  is  the  transmitter-to -target  range. 
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C.  RECENT  DEVELOPMENTS  IN  BISTATIC  AND  MULTISTATIC  RADAR 

Despite  research  progress  on  improving  the  detection,  classification,  and  location 
performance  of  passive  bistatic  radar  systems,  their  reliance  on  transmissions  of 
opportunity  and  the  restricted  geometry  has  limited  their  application.  This  fundamental 
requirement  continues  to  stimulate  research  and  experimentation  on  opportunistic 
illuminating  sources,  their  optimum  configuration,  applicability,  and  performance  in 
various  operational  scenarios.  As  with  all  radar  systems,  improvement  in  interference  and 
clutter  rejection,  target  detection,  classification,  and  location  and  tracking  accuracy  are 
current  areas  of  interest  in  the  field  of  passive  bistatic  radar  and  passive  multistatic  radar 
(PMR).  The  recent  research  topics  in  bistatic  and  multistatic  radar  can  be  classified  as 
bistatic  and  multistatic  system  configuration,  forward-scattering  radar,  and  multiple-input 
multiple-output  (MIMO)  radar. 

When  the  concept  of  bistatic  and  multistatic  radar  was  first  introduced,  dedicated 
radar  transmitters  were  used  as  transmission  sources  [11]  before  transmitters  of 
opportunity  were  employed  as  illuminating  sources.  Transmitters  of  opportunity  in  the 
very  high  frequency  (VHF)  and  ultra-high  frequency  (UHF)  band,  such  as  FM  broadcast, 
TV  broadcast,  DAB,  DVB-T,  and  cellular  network  signals  continue  to  be  common 
illuminating  sources  used  in  detecting  airborne,  land,  and  maritime  targets  [11-19]. 
Conversely,  studies  on  the  use  of  high  frequency  (HF)  band  signals  as  opportunistic 
transmitters  only  started  recently.  HF  signals  present  advantages  of  long  range  detection 
and  coverage,  propagation  beyond  the  radar  horizon,  and  improved  detection  of  stealth 
targets,  which  enhances  the  PBR/PMR’s  early  warning  capability  [20,  21].  A  list  of 
common  transmission  sources  and  their  typical  parameters  are  given  in  Table  2. 

Apart  from  terrestrial  sources  of  transmission,  there  has  been  recent  interest  in 
using  satellite  transmissions  to  detect  airborne  targets.  The  SABER-DEMO  platform 
recently  demonstrated  its  ability  to  detect  aircraft  passively  using  signal  processing 
techniques  to  process  weak  satellite  sources  of  transmission  [22], 
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Table  2.  Signal  parameters  for  typical  passive  radar  illumination  sources. 

From  [23], 


Transmission 

Frequency 

Modulation, 

Bandwidth 

P<Gt 

Power  Density  (Wm  2) 

<D-W 

4^Tj2 

HF  broadcast 

10-30  MHz* 

DSB  AM,  9  kHz 

50  MW 

-67  to  -53  dBW  m-2  at  n 
=  1000  km 

VHF  FM 
(analogue) 

-100  MHz 

FM,  50  kHz 

250  kW 

-57  dBW  m-2  at  n  =  100 
km 

UHF  TV 
(analogue) 

-550  MHz 

Vestigial-Sideband  AM 
(vision);  FM(sound), 

5.5  MHz 

1  MW 

-51  dBW  m'2  at  q  =  100 
km 

Digital  audio 
broadcast 

-220  MHz 

digital,  OFDM  220 
kHz 

10  kW 

-71  dBW  m'2  at  n  =  100 
km 

Digital  TV 

-750  MHz 

digital,  6  MHz 

8  kW 

-72  dBW  nU  at  n  =  100 
km 

Cellphone  base 
station  (GSM) 

900  MHz, 

1.8  GHz 

GMSK, 

FDM/TDMA/FDD  200 
kHz 

100  W 

-81  dBW  m'2  at  q  —  10 
0km 

Cellphone  base 
station  (3G) 

2  GHz 

CDMA  5  MHz 

100  W 

-81  dBW  m'2  at  =  100 
km 

*Appropriate  frequency  will  depend  on  time  of  day. 


Another  topic  of  interest  is  the  use  of  inverse  synthetic  aperture  radar  (ISAR) 
processing  techniques  with  forward-scattering  radar  (FSR)  for  target  detection  and 
parameter  extraction.  FSR  is  the  earliest  form  of  bistatic  radar  where  target  detection 
occurs  at  the  transmitter  to  receiver  baseline.  The  target’s  radar  cross  section  is  enhanced 
in  the  forward-scattering  configuration  due  to  Babinet’s  principle  [2],  Using  the  constant 
false  alarm  rate  (CFAR)  approach,  target  detection  and  parameter  extraction  can  be 
accomplished  in  real  time  [24],  An  experiment  conducted  at  Forsvarets  forskningsintitutt 
(FFI)  explores  the  difference  in  ISAR  ship  signatures  in  the  forward  and  back-scattering 
configurations.  The  results  from  the  study  show  that  forward-scattering  returns  produce 
more  accurate  ISAR  signatures  than  back-scattering  returns  as  the  RCS  for  forward¬ 
scattering  is  usually  stronger  than  back-scattering.  The  difference  between  forward  and 
back-scattering  ISAR  signatures  can  be  fused  to  improve  ship  identification  and 
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classification  [25],  Apart  from  studies  on  FSR  processing  techniques  to  improve  target 
detection  and  extraction,  considerable  research  has  been  conducted  on  ultra-wideband 
(UWB)  FSR  for  its  ability  to  reduce  sea  clutter  in  maritime  applications  [26,  27]. 

Recent  advances  and  interest  in  MIMO  radar  systems  can  be  attributed  to  its 
potential  for  detection  and  location  of  targets  in  bistatic  or  multistatic  configurations. 
Spatial  diversity  in  MIMO  radar  systems  exploits  the  differences  in  target  cross  section  in 
detecting  and  extracting  target  parameters  such  as  angle  of  arrival  (AO A)  and  Doppler 
frequencies  [28].  Furthermore,  coherent  processing  in  MIMO  systems  improves  target 
location  accuracy  [28,  29], 

As  with  all  studies  in  the  radar  domain,  current  bistatic  and  multistatic  research 
areas  are  motivated  by  the  need  to  improve  detection  capability  and  measurement 
accuracy  while  reducing  or  mitigating  the  effects  of  unwanted  interference. 

D.  THESIS  OBJECTIVE 

The  primary  objectives  of  this  thesis  are  as  follows: 

1.  Generate  a  MATLAB  model  that  computes  a  system’s  detection 
performance  given  the  target’s  path. 

2.  Examine  low-RCS  maritime  target  detection  coverage  and  performance 
using  multiple  pairs  of  bistatic  radar. 

3.  Investigate  low-RCS  maritime  target  location  accuracy  using  elliptical  and 
hyperbolic  target  location  methods. 

The  EM  simulation  software  FEKO  will  be  employed  to  model  and  analyze  the 
RCS  scattering  properties  of  the  maritime  target  while  MATLAB  will  be  used  to  model 
and  simulate  the  operating  environment  and  PBR  target  detection  and  location.  The 
MATLAB  model  is  validated  against  open  source  literature  by  using  simple  targets  with 
known  RCS  return  characteristics  prior  to  using  the  models  to  generate  simulation  results. 

The  MATLAB  model  generated  for  objective  (1)  requires  the  user  to  provide  an 
Excelfile  with  the  target’s  position  and  velocity  vector  components  at  each  time  step, 
after  which  detection  performance  parameters  will  be  computed.  This  model  will  be  a 
fully  automated  MATLAB  model  that  calls  FEKO  to  compute  precise  RCS  returns  given 
the  exact  incident  and  receive  angles.  Detection  coverage  results  from  objective  (2)  will 
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be  compared  against  findings  in  [4]  and  used  to  propose  the  optimal  PBR  configuration 
and  geometry  for  maximum  detection  coverage.  Target  location  estimation  results  from 
objective  (3)  will  be  used  to  examine  the  effects  of  PBR  geometry  on  location  error 
distribution. 

The  maritime  scenario  with  a  single  receiver  and  multiple  transmitters  depicted  in 
Figure  2  applies  to  all  models.  The  problem  setup  assumes  a  warship  (receiver)  equipped 
with  broadband  EW  receiver  and  direction  finding  (DF)  capability  deployed  to  monitor 
maritime  traffic  flow  in  the  Straits  of  Singapore.  The  targets  are  assumed  to  be  low-RCS 
targets  with  infrequent  transmissions  and/or  operating  with  LPI  radars,  while  civil  marine 
radars  (CMR)  on  merchant  ships  will  serve  as  opportunistic  illuminators  in  a  PBR 
configuration  [4],  Regulation  19  of  SOLAS  Chapter  V  requires  all  merchant  and 
warships  to  carry  automatic  identification  systems  (AISs),  which  share  information  on  the 
ship’s  identity,  position,  course,  speed,  navigation  status,  and  safety-related  information 
[30]. 


A  AxnA 
'A An. 

/  ./  "A  V  N 


Legend: 

Target  (Tgt) 


~50nmi 


- > 

- > 


Naval  Ship  (Rx) 

Merchant  Ships  (Tx) 

Maritime  AIS  and  Nav  Radar  Tx  (direct) 
Nav  Radar  Tx  (indirect) 


Figure  2.  Pictorial  representation  of  the  PBR  setup.  After  [4], 
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E.  THESIS  CHAPTER  OUTLINE 

Chapter  I  introduced  the  concept,  history,  and  recent  developments  in  passive 
bistatic  radar.  The  goals  and  end  products  of  the  thesis  are  also  detailed  here. 

Chapter  II  provides  the  theoretical  background  on  PBR.  The  corresponding 
parameters  employed  to  develop  the  necessary  MATLAB  model  are  introduced.  The 
bistatic  range  equations,  detection  contours,  and  mathematical  concepts  related  to  target 
location  and  uncertainty  are  covered. 

Chapter  III  presents  the  design  approach  in  modeling  the  scenario  and  the 
problem  setup  using  FEKO  and  MATLAB.  The  results  from  verifying  the  detection 
coverage  models  using  simple  targets  and  PBR  geometry  against  known  results  are 
covered  in  detail.  The  methodology  used  to  verify  target  location  model  is  also  discussed. 

Chapter  IV  uses  the  model  generated  to  examine  detection  coverage  for  different 
PBR  transmitter-target-receiver  geometries.  Simulation  results  are  compared  against 
findings  in  [4]  and  used  to  propose  PBR  configurations  and  geometry  for  optimal 
detection  coverage.  The  results  for  elliptical  and  hyperbolic  target  location  methods  and 
their  corresponding  uncertainty  ellipse  are  presented.  The  effects  of  transmitter-target- 
receiver  geometry  target  location  accuracy  are  also  explored. 

Chapter  V  summarizes  the  research  findings  and  suggests  further  work  to  improve 
the  models  and  multiple  PBR  detection,  parametric  extraction,  and  location  capability 
and  accuracy. 
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II.  BISTATIC  RADAR  THEORY 


A.  DEFINITION 

Bistatic  radar  refers  to  a  radar  system  where  the  transmitter  and  receiver  are  at 
sufficiently  different  locations  such  that  the  angles  or  ranges  from  those  locations  to  the 
target  are  significantly  different  [31].  The  basic  bistatic  configuration  and  parameters  are 
defined  in  Figure  1. 


B.  RANGE  EQUATION 


The  bistatic  radar  range  equation  (BRRE)  gives  the  received  power  at  Rx  as  a 
function  of  the  system  parameter,  target  scattering  properties,  and  engagement  geometry. 
Solving  the  BRRE  for  the  range  product  gives  [2] 


where 


K  =  (Rj-RflXmx 


PTGTGRA2aBF2F2Gf> 
(4*  Y  kTsBn(S  /  N)mnLTLR 


(1) 


RT  =  transmitter-to-target  range, 

Rr  =  receiver-to-target  range, 

PT  =  transmitter  power  output, 

Gt  =  transmitting  antenna  power  gain, 

Gr  =  receiving  antenna  power  gain, 

A  =  wavelength, 

crB  =  bistatic  target  cross  section, 

FT  =  pattern  propagation  factor  for  transmitter-to-target  path, 
Fr  =  pattern  propagation  factor  for  target-to-receiver  path, 

Gp  =  processing  gain, 
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k  =  Boltzmann’s  constant  (1.38x10  23  J/K), 

Ts  =  receiving  system  noise  temperature, 

Bn  =  noise  bandwidth  of  receiver’s  pre-detection  filter,  sufficient  to  pass  all 
spectral  components  of  the  transmitted  signal, 

(S  /  N)mm  =  signal-to-noise  power  ratio  required  for  detection, 

Lt  =  transmitting  system  losses  (>1)  not  included  in  other  parameters, 

Lr  =  receiving  system  losses  (>1)  not  included  in  other  parameters, 

K  =  bistatic  maximum  range  product. 

In  the  bistatic  range  equation,  the  maximum  range  product  RtRr  replaces  R2  in 
the  monostatic  range  equation  where  R  =  Rr  =Rr  is  the  monostatic  transmitter-to-target 

and  target-to-receiver  range.  The  difference  between  the  transmission  path  and  receiving 
path  results  in  significant  differences  between  monostatic  and  bistatic  radar  operation. 

One  of  the  differences  is  that  monostatic  contours  of  equal  signal  strength  are 
constant  range  circles,  while  detection  contours  for  bistatic  radar  are  defined  by  ovals  of 
Cassini.  An  oval  of  Cassini  is  defined  as  a  locus  of  points  where  the  product  of  the 
distance  from  two  fixed  points  is  constant.  Figure  3  shows  the  Cassini  oval  for  two  fixed 
points  ( Fj  and  F2 )  separated  by  a  distance  of  2c . 


Applying  the  concept  of  Cassini  ovals  to  the  bistatic  triangle  in  Figure  4  with 
baseline  L  and  range  product  RtRr  ,  an  expression  for  constant  signal-to-noise  (SNR) 
power  ratio  can  be  derived  by  writing  Eq.  1  as  [2] 


S/N  = 


K 

Wr 


(2) 


where  S  /  N  is  the  signal-to-noise  power  ratio  at  R,  and  Rr  ,  and  K  is  the  bistatic  radar 
constant 
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From  the  geometry  in  Figure  4,  RT  and  Rr  are  converted  to  polar  coordinates 


( r,0 )  using  the  law  of  cosines: 


Rj  —  (r2  +  )  +  vL,  cos  0, 

Rr  =  (^2  +  L/)  -  rL  cos  0, 


(4) 


4'  .  (5) 

where  the  origin  is  at  the  midpoint  of  the  baseline.  Substituting  Eq.  3,  Eq.4,  and  Eq.  5 
into  Eq.  2  gives  an  expression  that  defines  constant  SNR  contours  [2]: 


S/N  = 


_ K _ 

(r2+Z?/4)2  -r2L2  cos2  0 


(6) 


Signal-to-noise  ratio  contours  generated  using  Eq.  6  for  10  dB  <S/N<  30  dB  and 
K  =  30 L4  are  given  in  Figure  5. 


Figure  5.  Ovals  of  Cassini,  contours  of  constant  SNR  (dB),  with  K  =  30 L4 . 

After  [2], 
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Given  that  L  ,  ( Rr  +  Rr  )  and  0R  are  obtained  and  measured  from  the 


opportunistic  transmitter  and  receivers,  RT  and  Rr  are  computed  as  [2] 


R;  =  (R\  +L2  +  2RrL  sin  ffj'2, 

R  (Rr+RRf-L2 
2{Rt  +  Rr  +  L  sin  0 R ) 


Using  the  law  of  cosines  on  the  bistatic  triangle  in  Figure  4  yields 


P  =  cos 


-i 


fR2T+R2R 


r2\ 


2  RtRr 


(7) 

(8) 

(9) 


When  L  =  2 %[k  ,  the  oval  forms  a  lemniscate  with  cusp  at  the  origin.  The  ovals  of 
Cassini  in  Figure  5  define  three  operating  regions  for  bistatic  radar: 


1 .  L  >  2\fic  with  Rt»  Rr.  Receiver  centered  region. 

2.  L  >  2\fic  with  RK  «  Rr .  Transmitter  centered  region. 

3.  L  <  2 \[k  .  Cosite  region  or  receiver-transmitter-centered  region. 

In  cases  where  the  target  echo  signal  strength  is  weak,  non-coherent  pulse 
integration  performed  after  the  envelop  detector  increases  SNR  by  a  factor  of  \fN  where 
N  is  the  number  of  pulses  integrated.  Improvement  in  SNR  by  pulse  integration  is  a 
form  of  processing  gain.  The  number  of  pulses  integrated  over  a  period  of  t  seconds  is 
calculated  as 


where 


N -  PRF  xTOT xt 

^IdB  60 


TOT  = 


360  co. 


(10) 

(11) 


0MB  =  3  dB  azimuth  beamwidth, 


(oScan  =  scan  rate  (rpm), 

TOT  =  time  on  target, 

PRF  =  pulse  repetition  frequency, 

t  =  integration  period. 
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For  the  purpose  of  this  thesis,  the  minimum  difference  between  the  noise  floor  and  signal 
level  is  assumed  to  be  10  dB,  that  is,  minimum  SNR  is  10  dB. 


C.  TARGET  LOCATION  EQUATIONS 

1.  Bistatic  Radar  Trigonometry 

The  AOA  of  the  target  echo  signal  0R  and  target-to-receiver  range  Rr  are 

required  to  define  the  target’s  location  with  respect  to  the  receiver  in  a  PBR 
configuration.  The  AOA  of  the  echo  signal  can  be  measured  directly;  however,  the  target- 
to-receiver  range  cannot  be  measured  directly  and  needs  to  be  calculated  by  solving  the 
parameters  of  the  bistatic  triangle  (Figure  4). 

To  solve  for  target-to-receiver  range  Rr  and  the  rest  of  the  bistatic  triangle 
parameters  requires  measuring  and  knowledge  of  the  following: 

•  Baseline  range  from  transmitter  position(s)  and  receiver  position,  L  , 

•  AOA  of  target  echo  signal  at  the  receiver,  0R , 

•  Transmitter-to-target  and  target-to-receiver  range  sum,  (Rr  +  Rr  ) . 

The  range  sum  (Rr  +Rr)  can  be  estimated  using  the  direct  and  indirect  method  as 
illustrated  in  Figure  6.  In  the  direct  method,  the  receiver  measures  the  time  delay  A Trt 
between  the  reception  of  the  transmitted  pulse  and  the  target  echo.  The  range  sum  can 
then  be  expressed  as  a  function  of  the  time  delay  A Trt  and  the  baseline  range  [2] : 

(Rr  +Rr )  =  cATrt  +L.  (i2) 

In  the  indirect  method,  the  receiver  measures  the  time  delay  ATtt  between  the 
transmission  of  the  pulse  and  the  reception  of  the  target  echo.  The  range  sum  in  this  case 
is  a  function  of  the  time  delay  ATtt  [2] : 

(Rj.  +Rr)  =  cATtr  (13) 

The  indirect  method  requires  receiver  and  transmitter  clocks  to  be  synchronized  while  the 
direct  method  can  be  used  with  any  transmitter  configurations  given  LOS  between  the 
transmitter  and  receiver. 
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Referring  to  Figure  7,  the  direct  and  reflected  pulses  received  must  be  resolvable 
such  that 


^ reflected  ^  ^ direct  ^ 


(14) 


a)  Timing  Sequence 


TIME 


Pulse  emitted 
by  Tx 


Pulse  arrives  Pulse  arrives 

at  Rx  at  Tgt 


Tgt  echo 
arrives  at  Rx 


b)  Direct  Method 


4-  —  t'A  7^  4-  Z, 


AZ 


Start  Clock 


Stop  Clock 


c)  Indirect  Method 


R, 


+  Rlt-cAT,f 


AT, 


Start  Clock 


Stop  Clock 


Figure  6.  Timing  sequence  diagram  for  direct  and  indirect  method  for 
calculating  range  sum  ( Rr  +Rr)  .  From  [2], 
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Figure  7.  Reception  of  direct  and  reflected  pulses. 


To  derive  the  mathematical  relationship  between  the  bistatic  triangle  parameters, 
first  consider  the  elliptical  iso-range  contours  on  a  bistatic  plane  such  that  each  concentric 
ellipse  is  determined  by 

Rj.  +  Rr  2  a,  (15) 

where  a  is  the  semi-major  axis  length  of  the  ellipse  (Figure  3).  The  eccentricity  of  the 
ellipse  e  is  therefore  defined  as 


e  =  - 


2  a 
L 


e  =  ■ 


(Rr  +Rr) 


(16) 


(17) 


Given  the  measurement  of  (Rr+RR),  L ,  0R  and  using  the  law  of  cosines  on  the 


bistatic  triangle  in  Figure  4, 

R^=R2r+L2-  2  RrL  cos(90°  +  0R), 

R  (Rr+Rj^l 

R  2(RT +Rr+L  sin  dR)’ 

Rt  =  (R2r  +  L2  +  2  RrL  sin  0R  )1/2 . 

Substituting  the  range  sum  (Rr  +  Rr  )  in  Eqs.  19  and  20  using  Eq.  17  yields 


(18) 

(19) 

(20) 
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(21) 


Rr  = 


£(l~g2) 
2^(1  +  ^8016^) 


Rj  = 


L(e2  +  l  +  2esin  0R) 
2e(l  +  esin  0R) 


(22) 


Using  the  law  of  sines  on  the  bistatic  triangle  defines  the  relationship  between  the  range 
and  angle  values  as 


R„ 


Rr 


sin(90°  -0T)  sin(90°  +  6 R)  sin p 

Hence,  the  bistatic  angle  ft  is  expressed  as 


(23) 


P  =  sin 


P  -  sin 


f  t  h 


(24) 


(25) 


The  direct  and  indirect  method  of  measuring  range  sum  (RT  +Rr)  applies  to  all 

target  locations  except  when  the  target  is  in  a  forward-scattering  configuration  such  that 
the  target  lies  on  the  baseline  between  the  transmitter  and  receiver.  In  the  forward¬ 
scattering  PBR  configuration,  (Rt+Rr)  =  L  and  0R  =  -90  ,  making  Rr  in  Eq.  19 
indeterminate  [2], 


In  view  of  this  limitation  for  PBR  in  forward-scattering  configurations  and 
inaccuracies  arising  from  estimating  the  target’s  location  using  a  single  PBR  with 
erroneous  time  delay  and  AOA  measurements,  as  well  as  inaccurate  transmitter  and 
receiver  position  information,  the  following  sub-section  introduces  the  least-squares 
solution  for  fusing  bistatic  triangle  parameters  from  all  PBR  pairs. 


2.  Least-Squares  Intersection  of  Lines 

In  a  realistic  scenario,  AOA  and  range  information  for  each  PBR  pair  derived 
from  measurements  and  solving  trigonometric  equations  do  not  result  in  bearings 
intersecting  at  a  single  point  (Figure  8).  The  least-squares  solution  derived  in  [32]  finds  a 
point  that  minimizes  the  sum  of  perpendicular  distances  from  this  point  to  all  the  lines. 
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This  method  (i.e.,  the  one  that  solves  the  bistatic  triangle  parameters  and  estimates  the 
target  location  by  least  squares  bearing  intersection)  is  referred  to  as  the  elliptical  target 
location  method. 


Legend: 

X  Target 

-  Txl  bearing  from  Txl  position  and  0T  - Rx  bearing  from  Rr  and  0R  for  Txl 

-  Tx2  bearing  from  Tx2  position  and  0T  - Rx  bearing  from  Rr  and  0R  for  Tx2 


Tx3  bearing  from  Tx3  position  and  0T 


Rx  bearing  from  Rr  and  0R  for  Tx3 


Figure  8.  Least-squares  intersection  of  lines  solution  to  three  PBR  case. 


A  two-dimensional  line  is  described  by  a  point  on  the  line 


(26) 


and  its  corresponding  unit  direction  vector 


b  = 


brb  =  1 . 


(27) 
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The  squared  perpendicular  distance  from  a  point  p  to  a  line  as  illustrated  in  Figure  9  is 
expressed  as 


D( p;  a,  b)  =  ||(a  -  p)  -  ((a  -  p)rb)b||2 
=  (a-pf(I-bbr)(a-p), 
and  the  sum  of  squared  distance  for  K  lines  is 


(28) 


£>(p;a,b)  =  £z>(p;ay,by) 


7=1 


=Z( 


a.-p)r(I-b,b;)(a7-p). 


(29) 


Figure  9.  Perpendicular  distance  from  a  point  to  a  line.  From  [32], 


The  corresponding  objective  function  that  finds  the  “best-fit”  intersection  point  by 
minimizing  the  sum  of  squared  distances  for  all  lines  is 

argmin 


P  = 


Z)(p;A,B). 


(30) 


Taking  the  derivative  of  the  cost  function  with  respect  to  p 


dD 

dp 


X-2(I-bXXa,-p)  =  0, 


j= i 


Sp  =  q, 


21 


(31) 

(32) 


where 


S  =  Z(I-bX).  q  =  E(I-b,b>r 

j=1  j=l  (33) 

Solving  for  p  in  Eq.  32  gives  the  “best-fit”  point  of  intersection 

P  =  S"'l  (34) 

For  each  PBR  pair,  two  lines  are  defined  after  solving  for  their  bistatic  triangle 
parameters:  the  transmitter-to-target  bearing  and  the  target-to-receiver  bearing.  The 
transmitter-to-target  bearing  is  defined  by  the  transmitter  position  and  bistatic  triangle 
parameter^  while  the  target-to-receiver  bearing  is  defined  by  the  point  determined  by 

( 0r,Rr )  and  6 R .  In  a  forward-scattering  configuration,  where  Rr  is  indeterminate,  the 
target-to-receiver  bearing  is  defined  by  the  receiver  position  and  the  bistatic  triangle 
parameter  0R . 


3.  Hyperbolic  Target  Location 

Apart  from  the  elliptical  target  location  method  covered  in  Section  II.C.  1  and  its 
extension  in  Section  II.C.2,  a  target’s  position  can  also  be  estimated  using  a  hyperbolic 
location  technique.  To  derive  an  iterative  least-squares  method  of  estimating  the  target’s 
location  given  time  delay  measurement  from  the  direct  method  illustrated  in  Figure  6, 
consider  a  two-dimensional  PBR  pair  in  a  multi -bistatic  radar  scenario  (Figure  10). 
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Tgt 


Figure  10.  Single  PBR  in  multi-bistatic  radar  scenario. 

Given  the  coordinates  of  the  target  (xe,ye)  ,  receiver  position  (x 
transmitter  positions  (xtxi ,  yui ) ,  the  noiseless  time  delay  measurement  in  Eq. 


rewritten  as 

cAT  Rfi  +  Rr  L. 

hi  (x,  >  X,x,i  ’x^)  =  <J(xe  -  xtx,i  f + (ye  -  ytx,i  f 
+ \j(Xe  ~  Xrx  )2  +  (ye  ~yJ2 
-^rx-^f+iyrx-ytxj)2 

where 

i  = 

transmitter  number  i  =  1, 2,  •  •  • ,  K , 

K 

number  of  transmitters, 

R T,i  ~ 

Ith  transmitter-to-target  range, 

Rr  = 

target-to-receiver  range, 

L 

baseline  range, 

rx>yJ  and 

12  can  be 

(35) 


(36) 
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AT’  =  time  delay  between  zth  transmitter’s  direct  and  indirect  signal, 

Xe  =  target’s  position  vector,  [xt,  ye]T  , 


X 


ith  transmitter’s  position  vector,  [xtt  ;  yaJ  Jr , 
receiver  position  vector,  [xrx  yrv  ]r . 


When  noisy  time  delay  measurements  are  used  to  estimate  the  target  location,  Eq. 
35  is  written  as 

yi=c\Ti+ni 

=  *,(x„x„,x„)+«, 

where 


(37) 


y, 

A  Tt 

n: 


ith  noisy  time  delay  measurement, 


zth  noiseless  time  delay  measurement, 


=  zth  time  delay  measurement  error. 

Since  the  function  in  Eq.  36  is  a  non-linear  function  of  the  target,  receiver,  and 
transmitter  positions,  the  function  hi{Xe,Xtxi,\rx^j  will  be  linearized  by  a  Taylor  series 

expansion  about  an  initial  estimate  of  the  target’s  location  (xe0,ye0) .  By  retaining  the  first 
order  terms,  Eq.  37  can  be  written  as 


>•,-*,(  x, 0 ,  x„ ,  X„  )  =  Ax,  +  ft  Ay,  +  n„ 

dxe  8ye 


(38) 


where 


Axe 

= 

~Xe~ 

+ 

~  Xe0  ~ 

Je_ 

Je0_ 

(39) 

For  K  time  delay  measurements  ( K  PBR  pairs),  Eq.  38  can  be  represented  by  a 
linear  model 


Y  =  H  X  +  N  , 

(Xxl)  (Kx 2)  (2x1)  (Xxl) 
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(40) 


where 


Y  = 


'  y,-h,(X,0,XBJ,X„)  ' 

y2~K  (x,o.x«.x„) 

y  k  —  hK  (xg0 ,  x^ ,  x„  j 


,  H  = 


dl\ 

dhx 

dxe 

¥ 

dh2 

dh2 

dxe 

¥ 

dhK 

dhK 

dxe 

dye 

,  x  = 


Axe 

A ye 


(41) 


The  least-squares  solution  X  that  minimizes  the  sum  of  squares  of  difference  between  the 
measurements  and  the  estimated  function  is  defined  as 


X  =  (HrH)  HrY, 


(42) 


where 


X  = 


AXe 

at. 


The  estimated  target  location  in  the  current  iteration  is  therefore 


X,=xe0+x, 


(43) 


(44) 


such  that  the  target  location  estimate  in  the  current  iteration  is  used  as  the  initial  estimate 
Xe0  in  the  subsequent  iteration. 

D.  MEASUREMENT  AND  LOCATION  ERRORS 

The  theoretical  root  mean  square  (rms)  error  8M  of  a  radar  measurement  M  can 
be  expressed  as  [33] 

kM  kM 


dM  = 


yl2E/N0  4s7n 


(45) 


where  A:  is  a  constant  whose  value  is  of  the  order  of  one,  E  is  the  received  signal  energy, 
and  N0  is  the  noise  power  per  unit  bandwidth. 
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For  time-delay  (range)  measurements,  k  depends  on  the  shape  of  the  frequency 
spectrum  S(f) ,  and  M  is  the  rise  time  of  the  pulse.  For  angle  measurements,  &  depends 
on  the  shape  of  the  aperture  illumination  A(x) ,  and  M  is  the  beamwidth. 


1.  Time  Delay  (Range)  Measurements 

The  theoretical  rms  error  in  time  delay  measurements  dTR  for  a  rectangular  pulse 

with  pulsewidth  r  and  limited  by  bandwidth  B  is  approximately  [34] 

f  _  V/2 


dT„ 


V4  BE/N0, 


and  can  be  expressed  as  a  function  of  SNR: 


dT„ 


1 


b4sTn' 

This  assumes  Br  »1,  which  is  not  always  satisfied, 
rectangular  pulse  ( Bz  ~ 1 )  for  which  [33] 

r 


dTD 


2.1\fs  1  ~N 


(46) 


(47) 

A  more  accurate  model  is  a  quasi- 


(48) 


2.  Angle  Measurements 


The  theoretical  rms  error  of  AOA  measurements  dOR  for  an  antenna  with  uniform 


(rectangular)  amplitude  illumination  across  the  aperture  is  [33] 


ddR=- 


x/31 


0.628^„ 


7tD(2E/N0 )  (2  EIN0) 

and  can  be  expressed  as  a  function  of  SNR: 

=  O.62S0b 


1/2  9 


Ob  = 


4sTn  ’ 
0.88/1 
D 


where  D  is  the  antenna  dimension  and  0B  is  the  half-power  beamwidth. 


(49) 


(50) 

(51) 
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3.  Transmitter  and  Receiver  Position  Accuracy 

Transmitter  positions  on  merchant  ships  are  updated  by  differential  GPS  (DGPS) 
systems  and  made  available  to  surrounding  vessels  by  onboard  AIS  units.  There  are  two 
types  of  AIS  transceivers  [4]: 

1 .  Class  A  onboard  commercial  vessels 

2.  Class  B,  used  by  leisure  and  smaller  crafts 

For  the  purpose  of  this  thesis,  merchant  ships  are  assumed  to  be  equipped  with  Class  A 
AIS  that  broadcasts  the  vessel’s  unique  identification,  position,  course,  and  speed 
information  every  2  to  10  seconds  while  underway,  and  every  three  minutes  while  at 
anchor  at  a  power  level  of  12.5  W  [30], 

Receiver  positions  on  warships  are  also  provided  by  onboard  DGPS  systems.  In 
this  thesis,  the  U.S.  Coast  Guard’s  (USCG’s)  DGPS  service  accuracy  of  2  distance-root- 
mean-square  (drms)  [35]  is  used  to  model  transmitter  and  receiver  position  accuracies. 


4.  Receiver-to-Target  Range  Error 


As  explained  in  Section  C.l,  the  receiver-to-target  range  Rr  is  calculated  from  the 
range  sum  (Rt+Rr),  receiver  look  angle  6 R ,  and  baseline  L  measurements.  Assume 


that  the  measurement  errors  associated  with  range  sum  (RT  +Rr),  receiver  look  angle  0R , 

and  baseline  L  are  uncorrelated,  zero-mean  random  processes  having  Gaussian 
distribution  with  standard  deviation  equal  to  measurement  rms  error.  The  geometry 
dependent  root-sum-squared  (rss)  error  of  Rr  corresponding  to  Eq.  19  is  expressed  as  [2] 


d Rr  = 


OR,., 


“|2 


d  (Rt+Rr) 


d  (Rt+Rr) 


+ 


dR , 

dL 


d L 


+ 


dR, 


dOD 


(52) 


where  d(RT  +  Rr)  ,  dL  and  dOR  are  the  rms  errors  for  (Rt+Rr)  ,  L  ,  and  0R  , 
respectively,  and  determined  by  Eq.  45. 


From  Eq.  52,  the  components  of  the  rss  error  estimate  of  Rr  can  be  expressed  as  a 
function  of  eccentricity  e  and  receiver  look  angle  0R  [2]: 
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8Rr  1  +  6  +  2<?  sin  0 R 
8(Rt+Rr)~  2(l  +  esin6>fi)2  ’ 

8Rr  ( c  + 1)  sin  0R  +  2  c 


(53) 


dL  2(l  +  esin#R)2 
8Rr  _  L{ l-e2)  cos  0R 


(54) 


86r  2(1 +  e sin 6Rf  ’ 


(55) 


where  the  elliptical  iso-range  contour  eccentricity  is  e-L/(RT+RR )  and  each  partial 

derivative  component  defines  the  slope  of  the  error  surface  with  respect  to  each 
measurement  variable. 

E.  ERROR  ELLIPSE  PARAMETERS 

The  error  ellipse  provides  a  graphical  means  of  viewing  uncertainty  associated 
with  position  estimates.  The  error  ellipse  is  described  by  three  parameters:  (1)  semi-major 
axis,  a ;  (2)  semi-minor  axis,  b ;  and  (3)  orientation,  0  (Figure  1 1). 


v 


Semi-major  axis 


x 


Figure  1 1 .  Error  ellipse  parameters. 
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The  error  ellipse  parameters  are  described  by  the  covariance  matrix  of  the  target’s 
position  estimate,  Cx  y  ,  given  by 


= 


Pyexe  ^ ye  ^ xe 


Pxeye  ^ xe  ^ ye 

_ 2 


where  cr2  ,  cr2  and  are  the  variance  and  correlation  associated  with  the  target’s 
location  estimate. 

For  the  hyperbolic  target  location  method,  the  target  location  estimate  covariance 
matrix  is  expressed  as 

c«=[Jf  (57) 

where  C  is  the  covariance  matrix  for  time  delay  measurement  errors  and  J  „  is  the 
Jacobian  matrix  of  partial  derivations  in  Eq.  57: 

\an  0  0  0  ] 

n\ 

c  _  0  ^  0  0 

"  o  o’-,  o’ 


)  0 

d\ 

dhx 

1  ^ 

1  <8> 

dye 

dh2 

dh2 

dxe 

8ye 

dhK 

dhK 

dxe 

zye 

The  semi -major  axis  a  and  semi-minor  axis  b  are  expressed  as  a  function  of 
elements  in  Cx  y  as  [36] 


“^2 


(cr2  +cr2  )  +  .  —(cr2  +cr2  )  +(p  cr  cr  ) 

V  xe  ye)  \  a  \  xe  yd  \  ye  / 


b = + o  -  ■ 
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which  can  be  further  reduced  to  a  general  form 


a  =  4zlA 

b  =  s[zlA  \ 


(62) 


where  is  the  critical  value  of  the  chi-square  distribution  with  2  degrees  of  freedom 


and  evaluated  at  a  (Figure  12). 
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Figure  12.  Chi-square  probability  density  function  with  2  degrees  of  freedom. 
The  area  to  the  right  of  x2  critical  value  is  a  . 


There  are  two  forms  of  equations  used  to  compute  the  orientation  of  the  error 
ellipse  [36,  37]: 


6  =  —  tan  1 
2 


b  2  p  a  a  ^ 

r' xeye  xe  yt 

2  2 

v  y 


1.  ,-i 


f  cr2  -<j2  ^ 

9  =  —  tan"1  — * - 

2  j 


(63) 


(64) 
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where  Eq.  63  and  Eq.  64  represent  clockwise  and  anticlockwise  rotation  of  the  error 
ellipse  to  achieve  statistical  independence  (Figure  13). 


Figure  13.  Error  ellipse  rotation  to  achieve  statistical  independence.  From  [38]. 


To  understand  the  relationship  of  Monte  Carlo  simulations  of  target  locations  to 
their  corresponding  error  ellipses,  see  the  derivation  of  Eq.  62  and  Eq.  63  in  Appendix  A. 
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III.  MATLAB  AND  FEKO  MODELING 


A.  PROBLEM  SETUP 

The  maritime  scenario  with  a  single  receiver  and  multiple  transmitters  depicted  in 
Figure  2  will  be  employed  in  all  models.  The  problem  setup  assumes  a  warship  (receiver) 
equipped  with  broadband  EW  receiver  and  DF  capability  deployed  to  monitor  maritime 
traffic  flow  at  the  Straits  of  Singapore.  The  targets  are  assumed  to  be  low-RCS  targets 
with  infrequent  transmissions  and/or  operating  with  LPI  radars  while  civil  marine  radars 
(CMRs)  on  merchant  ships  will  serve  as  opportunistic  illuminators  in  PBR  configuration 
[4],  Regulation  19  of  SOLAS  Chapter  V  requires  that  all  merchant  ships  and  warships 
carry  automatic  identification  systems  (AISs)  which  share  information  on  the  ships’ 
identity,  position,  course,  speed,  navigation  status,  and  safety-related  information  [30]. 

The  target  detection  process  starts  with  the  warship  receiving  AIS  transmissions 
from  surrounding  transmitters,  allowing  the  warship  to  form  a  tactical  picture  with 
knowledge  of  the  transmitter  positions  and  the  corresponding  baseline  ranges.  Onboard 
EW  and  DF  receivers  will  identify  merchant  ships,  which  will  allow  separation  of  direct 
and  indirect  (target  echo)  transmissions  from  the  same  CMR.  The  time  delay  between 
reception  of  the  direct  and  indirect  transmissions,  DF  of  target  echo  and  baseline  range  is 
used  to  compute  the  target-to-receiver  range  in  all  PBR  configurations  except  in  a 
forward-scattering  configuration.  An  extension  to  this  elliptical  target  location  method  is 
to  define  transmitter-to-target  and  target-to-receiver  bearings  for  each  PBR  and  estimate 
the  target  location  by  finding  a  “best-fit”  intersection  via  least-squares.  An  alternative  to 
the  abovementioned  target  location  method  is  the  hyperbolic  target  location  approach, 
which  uses  time  delay  measurements  and  knowledge  of  the  transmitters  and  receiver 
position  to  estimate  the  target’s  position  iteratively. 

In  this  thesis,  technical  specifications  from  commercially  available  systems  are 
used  to  model  transmitters,  EW  and  DF  receivers.  Technical  specifications  of  UK-based 
Kelvin  Hughes’  MantaDigital  Radar  (see  Appendix  B)  in  the  S-band  and  X-band  are  used 
to  model  the  opportunistic  transmitters. 
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For  the  purpose  of  this  thesis,  merchant  ships  are  assumed  to  be  equipped  with 
Class  A  AIS  that  broadcasts  the  vessel’s  unique  identification,  position,  course  and  speed 
information  every  2  to  10  seconds  while  underway,  and  every  3  minutes  while  at  anchor 
at  a  power  level  of  12.5  W  [30].  Receiver  positions  on  warships  are  also  provided  by 
onboard  DGPS  systems.  For  the  simulation,  the  USCG’s  DGPS  service  accuracy  of  2 
dims  [35]  is  used  to  model  transmitter  and  receiver  position  accuracies. 

Onboard  the  warship,  the  EW  receiver  would  measure  and  extract  transmission 
parameters,  after  which  emitter  classification  and  identification  takes  place.  Integrated 
with  the  EW  receiver  is  the  DF  receiver  which  determines  the  AOA  of  transmission 
signals.  The  EW  receiver  specifications  used  in  the  models  are  based  on  UK-based 
Teledyne  Defence  Limited’s  high  performance  QR026  receiver  (see  Appendix  C).  This 
wideband  EW  receiver  operates  from  0.5  GHz  to  18  GHz  and  has  a  sensitivity  of  -70 
dBm.  The  DF  receiver  specifications  are  based  on  South  Africa’s  Poynting  Innovations 
(Pty)  Ltd’s  DF-AA0062  antenna  (see  Appendix  D),  which  has  coverage  from  20  MHz  to 
6  GHz. 

B.  FEKO  MODEL 

The  EM  simulation  software  FEKO  was  employed  to  model  and  analyze  the  RCS 
scattering  properties  of  the  maritime  target.  A  generic  frigate  is  chosen  as  the  maritime 
target.  A  three-dimensional  model  of  the  frigate  was  created  in  Computer  Simulation 
Technology  (CST)  Microwave  Studio  (MWS)  in  [4]  and  imported  into  FEKO  where  free 
space  monostatic  and  bistatic  RCS  returns  at  S-Band  and  X-Band  can  be  generated. 
Figure  14  shows  the  FEKO  mesh  model  and  a  picture  of  the  actual  frigate  on  which  the 
model  is  based,  while  Figure  15  shows  the  model’s  coordinate  system.  All  surfaces  are 
assumed  to  be  perfect  conductors. 
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Figure  14.  Frigate  FEKO  model  (top)  and  actual  RSN  Formidable-class  frigate 

(bottom;  from  [39]).  Side-profile. 

Figure  16  shows  the  frigate  monostatic  RCS  and  bistatic  RCS  at  10°  incident 
angle  at  3.05  GHz.  The  difference  in  monostatic  and  bistatic  RCS  can  vary  up  to  70  dB 
for  different  look  angles.  The  frigate’s  bistatic  RCS  returns  vary  with  bistatic  angle  (3 
and  is  larger  than  the  monostatic  RCS  in  the  forward  scatter  direction.  The  PBR 
configuration  capitalizes  on  this  characteristic  to  detect  low-RCS  targets  using 
opportunistic  transmissions. 
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Figure  15.  Frigate  FEKO  model  and  coordinate  system. 
Frigate  Monostatic  RCS  Frigate  Bistatic  RCS 


o 


Figure  16.  Frigate’s  monostatic  RCS  (dBsm)  (left)  and  bistatic  RCS  (dBsm) 
with  incident  angle  of  10°  (right)  at  3.05  GHz. 

Figure  17  shows  the  frigate’s  bistatic  RCS  returns  at  1°  and  0.1°  resolution, 

respectively,  against  the  azimuth  return  angle.  From  the  plots,  it  is  evident  that  small 

variations  in  azimuth  return  angles  may  result  in  significant  difference  in  RCS  returns. 

Simulations  using  the  same  model  in  FEKO  also  suggest  that  small  variations  in  the 

azimuth  incident  angle  will  result  in  significant  variations  in  RCS  returns.  Therefore,  to 

36 


ensure  that  accuracy  of  the  model  is  maintained,  it  would  be  ideal  to  call  FEKO  to 
compute  the  exact  RCS  returns  given  the  required  incident  and  receive  angles. 


Frigate  Bistatic  RCS 


Figure  17.  Frigate  bistatic  RCS  (dBsm)  with  incident  angle  of  10°  at  3.05  GHz. 

RCS  at  1°  resolution  (left)  and  0.1°  resolution  (right). 

C.  MATLAB  MODEL 

The  maritime  scenario  is  modeled  in  MATLAB  with  the  frigate  target  RCS  from 
FEKO.  Detection  coverage  and  target  location  methods  are  simulated  using  the 
MATLAB  models,  after  which  the  results  are  examined  to  propose  optimal  PBR 
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configurations  and  geometry,  as  well  as  to  investigate  their  effects  on  estimated  target 
location  uncertainty. 

1.  Detection  Coverage  Model  1 

The  original  MATLAB  detection  coverage  model  ( Detection  Coverage  Model  1 ) 
computes  the  received  SNR  given  the  target’s  position  information.  This  model  requires 
the  following  input  file  and  parameters  to  generate  a  detection  coverage  summary  plot: 

•  Target  path  information  in  an  Excelfile 

•  Target  CADFEKO  model  file  (.cfx  file) 

•  Target  default  EDITFEKO  file  (.pre  file) 

•  Transmitter  parameters 

•  Receiver  parameters 

The  target’s  path  information  Excelsheet  is  of  the  form  given  in  Table  3,  where 
the  model  is  required  to  execute  a  FEKO  run  for  each  PBR  pair  at  each  time  step.  The 
target’s  position  and  velocity  vector  at  each  time  step  can  represent  the  actual  target’s 
course  or  it  could  be  a  grid  of  equally  spaced  arbitrary  target  positions.  The  detection 
coverage  model  generates  the  received  SNR  along  the  target’s  path  for  the  former,  and  a 
detection  coverage  plot  (receiver  SNR  are  grid  points)  for  the  latter. 


Table  3.  Target  path  information  format  in  Excelsheet. 


Time 

Position  Vector 

Velocity  Vector 

*i 

fl 

zi 

V» 

Vzl 

t2 

*2 

y2 

Z2 

Vy2 

Vz2 

K 

Xn 

y« 

Zn 

V,n 

Kn 

Vz„ 

The  target’s  CADFEKO  and  EDITFEKO  files  are  required  for  FEKO  to  model 
and  generate  the  RCS  returns  from  the  target  for  given  incident  and  receive  angles.  The 
CADFEKO  file  contains  the  target’s  mesh  model  while  the  EDITFEKO  file  includes 
execution  parameters  such  as  frequency,  incident  angle,  and  receiver  angles.  For  each 
FEKO  run,  parameters  in  the  EDITFEKO  file  are  modified  prior  to  running  the  FEKO 
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Solver  function  on  the  CADFEKO  model.  The  RCS  returns  are  then  read  from  the  .out 
file.  The  transmitter,  EW,  and  DF  receiver  S-band  and  X-band  parameters  are  based  on 
technical  specifications  of  commercially  available  systems. 

Figure  18  gives  an  overview  of  Detection  Coverage  Model  1  for  each  transmitter- 
target-receiver  pair  at  each  time  step.  The  entire  MATLAB  model  and  its  interface  with 
Exceland  FEKO  are  automated.  The  model  begins  by  computing  the  incident  and  receive 
angles  at  the  target  given  the  transmitter’s,  receiver’s,  and  target’s  current  positions. 
Given  the  incident  and  receive  angles,  the  EDITFEKO  file  is  modified  to  indicate  the 
required  signal  incidence  and  receive  angles.  The  FEKO  solver  is  then  executed  on  the 
target’s  CADFEKO  model  to  generate  an  .out  file  containing  the  RCS  return  for  the 
current  transmitter-target-receiver  geometry. 


Figure  18.  MATLAB  Detection  Coverage  Model  1  flowchart. 


2.  Detection  Coverage  Model  2 

Figure  19  gives  an  overview  of  the  Detection  Coverage  Model  2  for  each 
transmitter-target-receiver  pair  at  each  time  step.  Unlike  Model  1,  Model  2  extracts  RCS 

returns  via  interpolation  given  bistatic  incident  and  receive  angle  from  a  pre-computed 
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RCS  table  of  the  form  given  in  Table  4.  This  reduces  the  simulation  duration  greatly  as 
the  model  does  not  need  to  execute  FEKO  for  each  angle. 


Figure  19.  MATLAB  Detection  Coverage  Model  2  flowchart 


The  Detection  Coverage  Model  2  also  computes  the  received  SNR  given  the 
target’s  position.  This  model  requires  the  following  input  file  and  parameters  to  generate 
a  detection  coverage  summary  plot: 

•  Target  path  information  in  an  Excelfile 

•  Target  RCS  table  in  an  Excelfile 

•  Transmitter  parameters 

•  Receiver  parameters 
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Table  4.  RCS  table  format  in  Excelsheet. 


The  target’s  path  information  Excelsheet  and  RCS  table  would  be  of  the  form 
given  in  Table  3  and  Table  4,  respectively.  Like  Model  1,  the  target’s  position  and 
velocity  vector  at  each  time  step  can  either  represent  the  actual  target’s  course  or  it  could 
be  a  grid  of  equally  spaced  arbitrary  target  positions.  The  detection  coverage  model 
generates  the  receive  SNR  along  the  target’s  path  for  the  former,  and  a  detection 
coverage  plot  (receiver  SNR  are  grid  points)  for  the  latter.  The  table  of  RCS  values  for 
various  incident  and  receiver  bistatic  angles  is  pre-computed  in  FEKO  using  the  target’s 
CADFEKO  model.  The  bistatic  angle  resolution  for  the  RCS  table  can  be  improved  by 
generating  a  higher  resolution  RCS  table  prior  to  simulating  detection  coverage.  Lastly, 
transmitter,  EW,  and  DF  receiver  S-band  and  X-band  parameters  are  based  on  technical 
specifications  of  commercially  available  systems. 

3.  Simulation  Duration  and  Accuracy  Trade-off 

Preliminary  simulations  using  Detection  Coverage  Model  1  presented  the  issue  of 
simulation  duration  and  accuracy  trade-off.  Referring  to  the  flowchart  of  Model  1  in 
Figure  18,  it  takes  approximately  4  seconds  for  FEKO  to  compute  the  RCS  each  time  the 
FEKO  Solver  is  called.  The  time  taken  is  acceptable  if  the  model  is  required  to  execute 
FEKO  via  the  model  only  a  few  times.  However,  in  generating  detection  coverage  plots, 
the  frigate  path  Excelfile  contains  equally  spaced  grid  points  across  an  area  of  interest.  In 


41 


this  case,  the  model  calls  FEKO  at  each  grid  point  for  each  PBR  pair.  For  example, 
consider  a  40  km  x  40  km  area  with  grid  points  spaced  1  km  apart  (Figure  20);  there  will 
be  41x41  =  1 681  grid  points.  If  there  is  one  opportunistic  transmitter  in  the  grid  area,  it 
will  take  approximately  1681  x  4  x  1=  6728  seconds  =1.9  hours  to  generate  the  detection 
coverage  for  the  grid  area.  However,  if  there  are  8  opportunistic  transmitters  in  the  grid 
area,  it  now  takes  approximately  1681  x  4  x  8  =  53,792  seconds  =  14.9  hours  to  generate 
the  detection  coverage. 
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Figure  20.  Detection  coverage  of  a  40  km  x  40  km  area  of  interest  grid  points  at 

1  km  resolution. 


The  large  number  of  simulation  cases,  compounded  with  the  need  to  simulate 
various  grid  sizes/resolutions  and  number  of  transmitters,  resulted  in  the  alternative 
detection  coverage  model  (Detection  Coverage  Model  2). 

Referring  back  to  Figure  17,  small  variations  in  return  angles  may  result  in 
significant  differences  in  RCS  return;  therefore,  Model  2  reduces  detection  coverage 
simulation  duration  at  the  expense  of  accuracy  due  to  low-RCS  resolution. 
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Figure  21  illustrates  the  simulation  duration  and  accuracy  trade  space  for  Model  1 
and  Model  2.  Assume  that  the  RCS  table  used  by  Model  2  has  bistatic  angle  resolution  of 
1  degree  and  its  current  position  in  the  simulation  trade  space  is  indicated  in  Figure  21. 
The  accuracy  of  Model  2  can  be  improved  with  minimal  increase  in  simulation  duration 
by  pre-computing  RCS  tables  with  higher  resolution  bistatic  angles.  Generating  a  high 
resolution  RCS  table  requires  longer  computing  duration  compared  to  a  lower  resolution 
RCS  table.  However,  the  additional  computing  duration  occurs  one  time  upfront  and 
greatly  reduces  the  simulation  duration  of  Detection  Coverage  Model  2. 

In  view  of  the  constraint  on  computing  speed,  Detection  Coverage  Model  1  and 
Detection  Coverage  Model  2  will  be  generated  but  Model  2  will  be  used  to  simulate 
target  detection  coverage  and  target  location  performance  for  the  rest  of  this  thesis. 


Figure  2 1 .  Simulation  duration  and  accuracy  trade  space  for  detection  coverage 
models. 
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4.  Target  Location  Model 

The  MATLAB  target  location  model  estimates  the  target’s  location  given  a  set  of 
erroneous  measurements  and  parameters.  This  model  requires  the  following  input  file  and 
parameters: 

•  Target  path  information  in  an  Excelfile 

•  Erroneous  time  delay,  DF  measurements 

•  Erroneous  transmitter  and  receiver  positions 

The  target’s  path  information  Excelsheet  would  be  of  the  form  given  in  Table  3. 
The  target’s  position  and  orientation  are  used  to  generate  measurements  and  compute  the 
SNR  at  the  receiver,  which  in  turn  is  used  to  model  rms  errors  in  the  measurements.  The 
theory  behind  modelling  measurement  errors  is  covered  in  Section  II.D.  For  the 
hyperbolic  target  location  method,  the  error  ellipse  associated  with  each  position  estimate 
is  computed. 

D.  MODEL  VERIFICATION 

Prior  to  using  the  Detection  Coverage  Models  to  generate  simulation  results,  it  is 
necessary  to  verify  the  model  against  bistatic  radar  principles  and/or  using  simple  targets. 
Here,  the  detection  coverage  models  were  verified  using  the  following  approaches: 

•  Target  with  constant  RCS 

•  Target  plane  at  various  orientations 

According  to  Eq.  2  and  Eq.  3  the  SNR  of  the  target  echo  depends  on  transmitter 
parameters,  receiver  parameters,  target  RCS  oB  ,  transmitter-to-target  range  R,  ,  and 

target-to-receiver  range  Rr  .  When  transmitter  parameters,  receiver  parameters,  and  target 
RCS  are  kept  constant,  the  SNR  depends  solely  on  RT  and  Rr  .  For  each  pair  of  constant 
Rt  and  Rr  ranges,  the  contour  of  constant  SNR  traces  out  an  oval  of  Cassini,  as  covered 
in  Section  II.B. 

To  verify  that  the  Detection  Coverage  Models  satisfy  bistatic  radar  principles, 
consider  a  target  whose  RCS  is  constant  regardless  of  aspect  angle  or  range.  The  resultant 
detection  coverage  should  exhibit  elliptical  SNR  contours,  as  depicted  in  Figure  5.  The 
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detection  coverage  and  contour  plot  generated  from  the  Detection  Coverage  Models  is 
shown  in  Figure  22.  The  results  agree  with  bistatic  radar  theory  such  that  the  SNR 
contour  values  satisfy  Eq.2. 


Detection  Coverage  SNR  Contours  (dB) 


Figure  22.  Detection  coverage  contours  for  constant  RCS  target. 

Apart  from  verifying  the  model  against  bistatic  radar  principles,  it  was  tested 
using  a  simple  plate  target  at  various  orientations.  The  FEKO  model  of  the  target  plate 
used  is  shown  in  Figure  23.  Figure  24  shows  the  S-band  RCS  for  the  target  plane  where 
the  angle  represents  the  angle  <f>  from  the  x-axis  horizontal  in  the  xy-plane  {6  =  90° )  of 
the  FEKO  model. 
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Figure  23.  Target  plate  FEKO  model. 


Target  Plane  RCS 

90 


Figure  24.  Target  plate  S-Band  azimuth  RCS  (dBsm)  at  cp=90°  incidence. 
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To  verify  that  the  MATLAB  Detection  Coverage  Model  is  sound  for  various 
target  orientations,  detection  coverage  plots  for  the  target  plate  are  generated  at  several 
rotations  of  the  plate  around  the  z-axis.  From  the  plate’s  RCS  in  Figure  24,  one  would 
expect  higher  RCS  and  hence  SNR  at  the  faces  of  the  plate,  that  is,  at  90°  and  270°.  In 
addition,  as  the  plate’s  RCS  returns  exhibit  symmetry  about  a  certain  axis,  the  detection 
coverage  plots  at  selected  target  orientations  demonstrate  symmetrical  characteristics.  See 
Appendix  E  for  the  target  plate  detection  coverage  plots  and  discussion  of  its  RCS. 

The  elliptical  and  hyperbolic  target  location  methods  outlined  in  Section  II.C  were 
verified  using  measurements  in  the  absence  of  errors.  Using  error-free  measurements, 
estimated  target  locations  have  zero  errors  for  both  target  location  methods. 

The  hyperbolic  target  location’s  associated  error  ellipse  semi-major  length,  semi¬ 
minor  length,  and  orientation  are  verified  against  Monte  Carlo  simulations.  For  example, 
it  was  verified  that  approximately  90%  of  the  Monte  Carlo  estimates  fall  within  a  90% 
confidence  ellipse. 
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IV.  SIMULATION  RESULTS 


In  this  chapter,  the  detection  coverage  for  various  PBR  transmitter-target-receiver 
geometries  is  examined  using  Detection  Coverage  Model  2,  and  PBR  configurations  and 
geometry  for  optimal  detection  coverage  are  proposed.  The  results  for  elliptical  and 
hyperbolic  target  location  methods  and  their  corresponding  uncertainty  ellipses  are 
presented.  This  chapter  begins  by  investigating  how  the  number  of  opportunistic 
transmitters,  their  ranges,  and  geometry  relative  to  the  target  and  receiver  affects 
detection  coverage  in  an  area  of  interest.  Detection  results  from  Detection  Coverage 
Model  1  and  Model  2  along  the  target’s  arbitrary  path  are  also  presented.  Target  location 
results  from  elliptical  and  hyperbolic  target  location  methods  are  then  presented  with 
their  corresponding  error  ellipse.  The  effects  of  transmitter-target-receiver  geometry  on 
each  target  location  method’s  spatial  error  distribution  are  explored.  All  simulation  plots 
and  results  are  included  in  Appendices  F  through  H. 

A.  DETECTION  COVERAGE 

1.  Bistatic  RCS 

To  better  understand  the  detection  coverage  plots  of  the  PBR  scenario  presented 
in  this  section,  consider  the  bistatic  RCS  of  the  target  (Figure  14)  with  incident  angle  10° 
at  3.05  GHz  and  9.41  GHz  in  Figure  25.  The  frigate’s  bistatic  RCS  returns  vary  with  the 
incident  and  receive  angles  with  peaks  occurring  at  certain  return  angles.  Overall,  they 
have  similar  return  patterns  in  S-band  and  X-band.  It  is  also  observed  that  RCS  will  have 
a  peak  in  the  forward-scattering  look  angle. 

Detection  coverage  plot  results  in  this  section  are  obtained  using  Detection 
Coverage  Model  2  with  a  RCS  table  of  1°  angle  resolution.  The  RCS  values  from  the  pre¬ 
computed  table  are  extracted  by  interpolation  if  incident  and/or  receive  angles  are  not 
integers.  Each  point  (x,y)  on  the  detection  coverage  plot  represents  the  total  received 

SNR  at  the  receiver  when  the  target  is  at  position  (x,y).  The  colour  at  each  point 
indicates  the  value  of  total  SNR  at  the  receiver. 
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Frigate  Bistatic  RCS  at  3.05  GHz  Frigate  Bistatic  RCS  at  9.41  GHz 

90  90 


Figure  25.  Target  bistatic  RCS  (dBsm)  with  incident  angle  10°  at  3.0  GHz  (left) 

and  9.41  GHz. 


2.  General  Observations 

Referring  to  the  detection  coverage  results  (see  Appendix  F),  the  SNR  of  a  target 
echo  at  the  receiver  strongly  depends  on  the  bistatic  RCS  characteristic  of  the  target. 
Higher  SNR  will  always  be  observed  when  the  target-to-receiver  LOS  lays  in  the  peak  of 
the  target’s  bistatic  RCS  returns.  Since  the  target’s  bistatic  RCS  exhibits  high  RCS 
returns  in  the  forward-scattering  look  angle  (Figure  25),  a  PBR  with  the  transmitter- 
target-receiver  in  the  forward-scattering  configuration  results  in  the  best  detection  at  the 
receiver  for  a  fixed  range.  It  is  also  evident  from  the  detection  coverage  results  that 
transmitter-target-receiver  in  the  back-scattering  configuration  would  result  in  good 
detection  with  target  echo  SNR  decreasing  as  the  target’s  range  to  receiver  increases. 
Figure  26  highlights  target  positions  in  the  forward  and  back-scattering  configuration  on 
a  detection  coverage  plot  where  the  transmitter  is  positioned  at  (-5000,0)  and  receiver 
at  (0,0) . 


50 


Detection  Coverage  SNR  (dB) 


3 


2 


1 


E 


fo 

■R 

a 


-2 


-3 


-4 

-4-3-2-101234 

Distance  (m)  x  ig* 


Figure  26.  Target  positions  in  forward-scattering  and  back-scattering 
configuration  on  detection  coverage  plot. 

From  the  detection  coverage  results  (see  Appendix  E),  it  was  observed  that  target 
echo  SNR  is  generally  high  when  the  target  is  within  the  “net”  of  transmitters-receiver 
configuration  with  the  highest  SNR  occurring  when  the  PBR  is  in  the  forward-scattering 
configuration.  This  observation  agrees  with  the  findings  in  [4]  which  state  that  the 
optimal  transmitters-receiver  constellation  is  to  have  transmitters  at  every  quadrant  with 
respect  to  the  receiver.  Outside  the  “net”  of  transmitters-receiver,  detection  is  usually 
better  when  the  PBR  is  in  the  back-scattering  configuration.  Across  all  PBR  geometries, 
detection  gaps — areas  with  the  low  target  echo  SNR  at  the  receiver — occur  when  the 
target  is  outside  the  transmitter-receiver  cluster  or  away  from  bearings  extending  from  the 
transmitter-receiver  baselines  (Figure  27). 
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Detection  Coverage  SNR  (dB) 


Figure  27.  Detection  gaps  and  bearings  extending  from  Tx-Rx  baseline  in 

detection  coverage  plot. 


3.  Vary  Number  of  Transmitters 

The  number  of  available  opportunistic  transmitters  affects  the  detection  coverage 
area  and  strength  of  the  target  echo  at  the  receiver.  Figure  28  shows  that  overall  detection 
coverage  improves  as  the  number  of  opportunistic  transmitters  increases.  The  result 
presented  in  this  thesis  assumes  that  target  echo  is  received  from  all  transmitters  shown  in 
the  detection  coverage  plot.  The  percentage  of  grid  points  with  received  SNR  greater  than 
10  dB  within  a  10  km  radius  from  the  receiver  for  a  different  number  of  transmitters 
around  the  receiver  is  given  in  Table  5  and  0.  For  a  target  orientation  of  0°,  detection 
coverage  increases  from  17.0%  to  61.6%  as  the  number  of  opportunistic  transmitters 
increase  from  one  to  eight  (Table  5). 
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Figure  28.  Detection  coverage  plots  at  S-Band  for  with  target  plate  at  0° 
orientation  as  the  number  of  transmitters  varies. 

The  detection  statistics  presented  in  Table  5  and  0  suggest  that  S-band 
opportunistic  transmitters  result  in  better  detection  than  X-band  opportunistic  transmitters. 
This  is  attributed  to  the  variation  in  RCS  returns  in  X-band  compared  to  S-band. 
Referring  to  Figure  25,  peak  RCS  at  X-Band  is  higher,  but  narrow  lobes  result  in  low 
probability  of  intercepting  a  lobe.  Nevertheless,  detection  coverage  patterns  are  consistent 
between  S-band  and  X-band  transmissions.  As  the  difference  in  detection  coverage  is  not 
significant  in  detection  coverage  plots,  subsequent  sections  of  this  thesis  are  focused  on 
results  from  S-band  transmissions. 
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Table  5.  Percentage  of  grid  points  with  SNR  greater  than  10  dB  within  a  10 
km  radius  centered  at  the  receiver.  Target  orientation  at  0°. 


Target  at  0°  orientation 

#  of  Tx 

%  SNR  >  10  dB  in  10  km  radius 

S-Band 

X-Band 

1 

17.0% 

15.1% 

2 

33.8% 

29.5% 

3 

42.0% 

36.4% 

4 

47.5% 

40.3% 

8 

61.6% 

54.8% 

Table  6.  Percentage  of  grid  points  with  SNR  greater  than  10  dB  within  a  10 
km  radius  centered  at  the  receiver.  Target  orientation  at  -90°. 


Target  at  -90°  orientation 

#  of  Tx 

%  SNR  >  10  dB  in  10  km  radius 

S-Band 

X-Band 

1 

16.4% 

15.7% 

2 

33.4% 

29.5% 

3 

40.0% 

35.4% 

4 

47.5% 

40.3% 

8 

61.6% 

54.8% 

4.  Vary  Transmitter  Range 

The  transmitter  to  receiver  range  determines  the  size  of  the  transmitter-receiver 
“net”  and  forward-scattering  area.  Detection  coverage  plots  in  Figure  29  suggest  that  the 
area  of  detection  coverage  increases  as  the  transmitter  to  receiver  ranges  increase; 
however,  the  highest  target  echo  SNR  in  the  grid  decreases  with  increase  in  range.  This 
observation  is  consistent  with  the  findings  in  [40]  with  regards  to  sensitivity  dependency 
on  transmitter  range  and  transmitter-receiver  geometry.  The  percentage  of  grid  points 
with  received  SNR  greater  than  10  dB  within  a  20  km  radius  from  the  receiver- 
transmitter  range  of  5  km,  10  km,  and  15  km  are  listed  in  Table  7.  For  target  orientation 
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of  0°,  detection  coverage  increases  as  the  transmitter  to  receiver  range  increases  from  5 
km  to  10  km.  However,  the  detection  coverage  decreases  when  the  range  further 
increases  to  15  km  unless  more  than  four  transmitters  are  present. 


Figure  29.  Detection  coverage  plots  at  S-Band  with  target  plate  at  0°  orientation 

as  the  transmitter  range  to  receiver  increases. 
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Table  7.  Percentage  of  grid  points  with  SNR  greater  than  10  dB  within  a  20 
km  radius  centered  at  the  receiver.  Target  orientation  at  0°. 


%  SNR  >  10  dB  in  20  km  radius 

. . ^ . T x-Rx  Range 

#ofTx 

5  km 

10  km 

15  km 

1 

7.3% 

8.7% 

5.9% 

2 

14.4% 

17.3% 

12.7% 

4 

19.7% 

28.8% 

19.9% 

8 

31.2% 

37.1% 

38.0% 

5.  Vary  Transmitter-Target-Receiver  Geometry 

Two  types  of  transmitter-target-receiver  geometries  were  explored:  (1)  receiver- 
centered  geometry,  as  seen  in  Figure  30,  and  (2)  transmitter-clustered  geometry,  as  seen 
in  Figure  31.  In  a  receiver-centered  geometry,  all  observations  made  earlier  in  this 
chapter  hold.  High  target  SNR  is  received  when  the  target  is  in  the  transmitter-receiver 
“net”  with  the  best  detection  occurring  when  the  transmitter-target-receiver  is  in  the 
forward-scattering  configuration.  Outside  the  transmitter-receiver  “net,”  detection  is  the 
best  when  the  transmitter-target-receiver  is  in  the  back-scattering  configuration.  Apart 
from  areas  that  satisfy  a  back-scattering  configuration,  it  is  observed  that  there  are 
additional  areas  that  exhibit  good  detection  coverage  as  a  result  of  the  target’s  bistatic 
RCS  patterns.  Multiple  PBR  capitalizes  on  this  to  allow  detection  of  low-RCS  targets. 

In  a  transmitter-clustered  geometry,  good  detection  is  experienced  when  the  target 
lies  in  the  cluster  of  transmitters  or  in  the  beam  extending  from  the  receiver  to  the  cluster 
of  transmitters.  Outside  the  transmitter-receiver  cluster,  detection  is  slightly  better  at 
areas  that  satisfy  back-scattering  configuration  and  “lines”  that  are  shaped  like  a 
hyperbolic  curve — hyperbolic  contours  are  typical  of  time  difference  of  arrival  (TDOA) 
measurements  [37],  However,  for  a  four  opportunistic  transmitter’s  scenario,  the  received 
target’s  echo  is  still  below  10  dB  outside  the  transmitter-receiver  cluster. 
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Figure  30.  Detection  coverage  at  S-Band  for  4  Tx  in  receiver-centered  geometry 

and  target  plate  at  30°  orientation. 
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Figure  3 1 .  Detection  coverage  at  S-Band  for  4  Tx  in  transmitter-clustered 
geometry  and  target  plate  at  30°  orientation. 
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Up  to  this  point,  all  transmitters  in  a  scenario  were  placed  in  a  predictable  fashion 
within  an  area  of  interest.  The  next  series  of  detection  coverage  plots  randomly  positions 
20  transmitters  at  a  range  of  5  to  20  km  from  the  receiver.  Earlier  observations  in  regard 
to  areas  of  good  detection  coverage  are  evident  in  the  detection  coverage  plots.  However, 
strips  of  detection  gaps — positions  with  the  lowest  target  echo  SNR  at  the  receiver — are 
evident  in  Figure  32,  Figure  33,  and  Figure  34. 
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Figure  32.  Detection  coverage  at  S-Band  for  20  Tx  randomly  positioned  around 
Rx  at  5-30  km  range  and  target  plate  at  0°  orientation. 

The  direction  and  location  of  the  detection  gaps  are  determined  by  the  target’s 
velocity  vector  and  hence  its  orientation.  Referring  to  Figure  32  where  the  target  velocity 
vector  is  horizontal  and  points  to  the  right  and,  hence,  the  target  orientation  is  0°  across 
all  grid  points,  detection  gaps  will  present  themselves  as  horizontal  strips.  Applying  this 
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explanation  for  target  orientations  of  -90°  and  45°,  detection  strips  would  be  vertical  and 
angled  at  45°,  respectively  (Figure  33  and  Figure  34).  The  strips  of  detection  gaps  are 
observed  to  be  present  when  the  target’s  heading  is  towards  the  transmitters-receiver 
cluster. 
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Figure  33.  Detection  coverage  at  S-Band  for  20  Tx  randomly  positioned  around 
Rx  at  5-30  km  range  and  target  plate  at  -90°  orientation. 
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Figure  34.  Detection  coverage  at  S-Band  for  20  Tx  randomly  positioned  around 
Rx  at  5-30  km  range  and  target  plate  at  45°  orientation. 

6.  Target  Path  Detection 

Detection  Coverage  Model  1  and  Model  2  compute  detection  performance  given 
the  target’s  path,  with  Model  1  using  more  accurate  RCS  values  from  FEKO  and  Model  2 
extracting  RCS  values  from  a  pre-computed  table  given  bistatic  incident  and  receive 
angles.  The  trade-off  between  accuracy  and  simulation  duration  illustrated  in  Section 
III.C.3  limits  the  use  of  Model  1  as  a  tool  to  study  detection  performance.  Nevertheless, 
target  path  detection  performance  results  from  both  models  are  included  in  Appendix  G. 

Five  arbitrary  target  paths  in  Figure  35  are  used  to  generate  target  path  detection 
performance  as  the  target  moves  through  the  area  of  interest.  The  received  target  echo 
SNR  as  the  target  transverses  through  the  area  of  interest  is  given  in  Figure  36.  The 
detection  performance  results  here  align  with  previous  observations,  where  detection 
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coverage  is  good  within  the  transmitters-receiver  cluster  and  in  the  back-scattering 
configuration  outside  the  transmitter-receiver  cluster. 
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Figure  35.  Five  target  paths  used  to  generate  target  path  detection  performance 

results. 


Figure  36.  Detection  coverage  along  target’s  path  at  S-Band  for  4  Tx  (left)  and 
20  Tx  (right).  RCS  extracted  from  pre-computed  RCS  table. 


61 


B.  TARGET  LOCATION  ESTIMATION 

1.  General  Observations 

To  examine  target  location  estimation  accuracy  for  elliptical  and  hyperbolic  target 
location  methods,  results  from  both  methods  are  presented  in  target  location  error  color 
plots  and  estimated  target  location  scatter  plots.  In  a  target  location  error  plot,  the  color  at 
each  point  (x,y)  represents  the  error  in  location  estimate  when  the  true  target  is  located 

at  (x,y).  For  estimated  target  location  scatter  plots,  the  estimated  target  location  for  each 

round  of  Monte  Carlo  simulation  is  marked  such  that  spatial  error  distribution  for  each 
location  method  can  be  investigated. 

Referring  to  the  target  location  error  color  plots  where  measurement  errors  are 
independent  of  SNR  (see  Appendix  H),  the  hyperbolic  location  method  generally 
estimates  target  location  more  accurately  compared  to  the  elliptical  location  method  as  it 
uses  fewer  measurement  sources.  The  hyperbolic  method  uses  time  delay  measurements 
and  transmitters  and  receiver  position  data,  while  the  elliptic  method  uses  DF 
measurements,  time  delay  measurements  and  transmitter  and  receiver  position  data. 
When  measurements  with  errors  independent  of  SNR  are  used,  it  was  observed  that  the 
effects  of  dilution  of  precision  on  target  location  estimates  are  less  pronounced  for  the 
hyperbolic  method  (Figure  37). 

To  better  model  errors  present  in  measurements  used  to  estimate  the  target’s 
location,  target  location  models  will  include  SNR-dependent  errors  in  the  DF  and  time 
delay  measurements,  as  outlined  in  Section  II.D.  Simulating  target  location  accuracy  with 
various  combination  of  erroneous  and  error-free  measurements  suggest  that  DF 
measurement  is  the  dominant  source  of  error  in  the  elliptical  method,  while  time  delay 
measurement  is  the  dominant  source  of  error  in  the  hyperbolic  method.  It  was  also 
observed  that  the  effect  of  2  meters  dims  error  in  transmitter  and  receiver  position  data 
are  negligible  compared  to  SNR-dependent  DF  and  time  delay  measurements. 
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Figure  37.  Dilution  of  precision  comparison  between  elliptical  and  hyperbolic 
methods  using  SNR-independent  measurement  errors. 

Referring  to  the  target  location  error  results  where  measurements  with  SNR- 
dependent  measurement  errors  are  used  (Appendix  H),  the  hyperbolic  method,  which 
uses  fewer  sources  of  measurement  and  data,  produces  more  accurate  location  estimates 
compared  to  the  elliptical  method.  However,  overall  location  accuracy  is  poor  due  to  low 
SNR  and  large  DF  and  time  delay  measurement  errors. 

The  results  from  the  target  location  accuracy  color  plots  mirror  the  observation 
made  in  detection  coverage  plots.  Target  location  accuracy  is  best  when  the  target  lies 
within  a  transmitters-receiver  cluster  in  receiver-centered  geometry  (Figure  38). 
Increasing  the  transmitter  to  receiver  range  effectively  increases  the  transmitters-receiver 
cluster  area  but  reduces  the  overall  location  accuracy  within  the  cluster. 
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Figure  38.  Target  location  errors  using  hyperbolic  method  (S-Band)  with  SNR- 

dependent  measurement  errors. 

2.  Error  Ellipse  of  Target  Position  Estimate 

To  examine  the  spatial  error  distribution  of  target  position  estimates,  100  Monte 
Carlo  simulations  are  performed  for  each  transmitters-target-receiver  geometry.  The 
target  position  estimate  from  each  Monte  Carlo  simulation  is  marked  on  a  two- 
dimensional  plot  to  produce  a  target  location  scatter  plot.  For  the  hyperbolic  target 
location  method,  parameters  of  each  estimate’s  error  ellipse  can  be  calculated  for  various 
levels  of  confidence.  The  equations  used  to  compute  the  error  ellipse  semi-major  axis 
length,  semi-minor  axis  length,  and  orientation  is  covered  briefly  in  Section  II.E  and 
derived  in  Appendix  A. 

Transmitter-receiver  geometries  can  be  broadly  classified  into  two  groups:  (1) 
target  within  Tx-Rx  cluster  and  (2)  target  outside  Tx-Rx  cluster. 
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Figure  39  and  Figure  40  show  the  target  location  scatter  plot  from  the  elliptical 
and  hyperbolic  methods  respectively  when  group  (1)  geometry  occurs.  The  estimates 
from  all  100  Monte  Carlo  simulations  cover  an  approximately  circular  area,  suggesting 
that  the  extent  of  target  estimate  uncertainty  in  the  x  and  y  direction  is  similar.  The  area 
covered  by  the  scatter  points  is  referred  to  as  the  error  ellipse  for  position  estimate.  The 
error  ellipse  of  estimates  from  the  hyperbolic  method  is  smaller  than  that  from  the 
elliptical  method.  This  suggests  that  the  hyperbolic  method  estimates  target  position  more 
accurately. 

Similarly,  Figure  41  and  Figure  42  show  the  target  location  scatter  plot  from  the 
elliptical  and  hyperbolic  methods  respectively  when  group  (2)  geometry  is  considered. 
The  estimates  from  all  100  Monte  Carlo  simulations  in  this  case  form  an  ellipse, 
exhibiting  DOP.  The  resultant  error  ellipse  from  using  the  elliptical  and  hyperbolic 
method  differs  greatly  with  their  semi-major  axis  appearing  approximately  orthogonal. 
To  understand  how  measurement  errors  translate  to  uncertainty  in  the  target  location 
estimates,  consider  the  graphical  illustration  in  Figure  44  through  Figure  47. 


Estimated  Target  Locations  using  Elliptical  Method 
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Figure  39.  Target  position  estimate  scatter  plot  from  elliptical  method  for  target 

in  Tx-Rx  cluster. 
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Estimated  Target  Locations  using  Hyperbolic  Method 
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Figure  40.  Target  position  estimate  scatter  plot  from  hyperbolic  method  for 

target  in  Tx-Rx  cluster. 
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Estimated  Target  Locations  using  Elliptical  Method 
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Figure  4 1 .  Target  position  estimate  scatter  plot  from  elliptical  method  for  target 

outside  Tx-Rx  cluster. 
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Estimated  Target  Locations  using  Hyperbolic  Method 
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Figure  42.  Target  position  estimate  scatter  plot  from  hyperbolic  method  for 

target  outside  Tx-Rx  cluster. 

In  the  elliptical  location  method,  bistatic  triangle  parameters  are  computed  using 
DF  and  time  delay  measurements  as  well  as  transmitters  and  receiver  position  data.  The 
transmitter-to-target  range  RT ,  receiver  to  target  range  Rr  ,  and  DF  measurements  are 

used  to  define  several  bearings  which  would  intersect  at  the  target’s  true  position  if  all 
measurements  are  error-free.  However,  in  reality,  these  bearings  do  not  intersect  at  a 
single  point,  hence  the  elliptical  method  finds  “best-fit”  point  of  intersection  and  returns 
this  point  as  the  estimated  target  position.  Graphically,  uncertainty  associated  with  each 
transmitter-to-target  or  receiver-to-target  bearing  is  defined  as  a  sector  extending  from 
the  transmitter  or  receiver  towards  the  target’s  position.  The  intersection  of  all  sectors 
approximates  the  uncertainty  associated  with  the  target  position.  For  geometry  (1),  the 
target  position  estimates  from  the  elliptical  method  results  in  a  relatively  circular  and 
small  error  ellipse  (Figure  43),  while  for  geometry  (2),  the  error  ellipse  is  eccentric  with 
the  semi-major  axis  aligned  with  the  receiver’s  LOS  (Figure  44). 
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Figure  44.  Uncertainty  area  associated  with  position  estimates  from  elliptical 

method  for  target  outside  Tx-Rx  cluster. 


68 


For  the  hyperbolic  target  location  method,  target  position  is  estimated  by 
minimizing  the  sum  of  squares  of  difference  between  time  delay  measurements  and  the 
estimated  function.  The  range  equivalent  time  delay  measurement  is  defined  by  Eq.  12 
and  rewritten  as 

ciJ„=(Rr+R,)-L.  (65) 

For  a  transmitter-target-receiver  pair,  baseline  L  is  constant;  therefore,  the  time 
delay  measurement  will  depend  on  the  range  sum  (RT  +  Rr)  ■  As  in  Eq.  15,  range  sum 

equals  twice  the  semi-major  axis  length.  From  the  definition  of  an  ellipse,  the  iso-contour 
of  the  constant  time  delay  measurement  is  elliptical  with  the  transmitter  and  receiver 
position  as  ellipse  foci  (Figure  45). 
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Figure  45.  Elliptical  iso-contours  of  constant  time  delay  measurements  with 

transmitter  and  receiver  at  ellipse  foci. 


In  the  hyperbolic  target  location  method,  an  erroneous  time  delay  measurement 
from  a  transmitter-target-receiver  pair  defines  an  ellipse  with  “thickness”  in  which  the 
“thickness”  is  determined  by  the  amount  of  error  in  the  time  delay  measurement.  The 
intersection  of  all  ellipses  approximates  the  uncertainty  associated  with  the  target  position. 
For  geometry  (1),  the  target  position  estimates  from  the  hyperbolic  method  result  in  a 
relatively  circular  and  small  error  ellipse  (Figure  46),  while  for  geometry  (2),  the  error 
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ellipse  is  eccentric  with  the  semi-major  axis  perpendicular  to  the  line  extending  from  the 
target  towards  the  receiver  (Figure  47). 


Figure  47.  Uncertainty  area  associated  with  position  estimates  from  hyperbolic 

method  for  target  outside  Tx-Rx  cluster. 
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Comparing  the  scatter  plots  from  the  elliptical  and  hyperbolic  methods  in  Figure 
41  and  Figure  42,  respectively,  it  is  evident  that  uncertainty  associated  with  position 
estimates  from  the  hyperbolic  method  is  smaller.  The  error  distribution  of  position 
estimates  from  elliptical  and  hyperbolic  methods  appear  to  be  orthogonal  and  are  more 
pronounced  in  geometries  that  exhibit  DOP;  that  is,  the  semi-major  axis  of  the  error 
ellipse  of  estimates  from  the  elliptical  method  is  approximately  aligned  with  the  semi¬ 
minor  axis  of  the  error  ellipse  of  estimates  from  the  hyperbolic  method,  and  vice  versa.  A 
possible  option  to  improve  target  position  estimates  is  to  combine  the  elliptical  and 
hyperbolic  method  in  such  a  way  that  capitalizes  on  precision  in  their  respective  semi¬ 
minor  axis  direction. 

Every  target  position  estimate  is  accompanied  by  a  measure  of  confidence.  The 
method  of  computing  the  error  ellipse  parameters  outlined  in  Section  II.E  and  Appendix 
A  is  an  indication  of  the  estimate’s  accuracy.  The  error  ellipse  shows  how  measurement 
errors  translate  to  the  estimate’s  uncertainty  given  certain  geometry  with  the  ellipse  semi¬ 
major  and  semi -minor  axes  being  directions  with  greatest  and  least  deviation.  The  size  of 
the  ellipse  is  determined  by  the  critical  chi-square  value  %l  such  that  the  area  to  the  left 
(l-a)xl00%  is  its  corresponding  confidence  level.  For  a  90%  confidence  error  ellipse, 
the  critical  chi-square  value  is  4.60,  and  it  will  enclose  approximately  90%  of  the  Monte 
Carlo  position  estimates.  For  the  transmitters-target-receiver  geometry  in  Figure  42,  the 
error  ellipse  at  various  confidence  levels  for  hyperbolic  position  estimates  is  shown  in 
Figure  48.  It  was  verified  that  the  number  of  Monte  Carlo  estimates  that  fall  within  the 
ellipse  agrees  with  the  confidence  level  of  the  ellipse. 
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Figure  48.  Hyperbolic  target  location  estimate’s  error  ellipse  at  50%,  70%,  90%, 

and  99%  confidence  levels. 
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V.  CONCLUSION 


A.  SUMMARY  OF  FINDINGS 

This  thesis  explored  two  areas  pertaining  to  the  use  of  multiple  passive  non- 
cooperative  opportunistic  transmitters  in  target  detection  and  location.  To  investigate 
these  PBR  issues,  detection  coverage  and  target  location  estimation  models  have  been 
developed  using  FEKO  and  MATLAB  software.  The  models  were  used  to  simulate 
detection  performance  and  location  estimation  accuracy  for  various  transmitter-target- 
receiver  configurations  in  a  maritime  scenario  with  a  low-RCS  target. 

The  proposed  PBR  approach  offers  potential  advantage  in  the  detection  of 
stealthy,  low-RCS  targets  which  are  designed  to  minimize  monostatic  radar  echoes.  RCS 
returns  of  stealthy  ships  can  be  sufficiently  large  at  certain  return  directions,  allowing 
PBR  to  capitalize  on  this  characteristic  to  detect  low-RCS  targets  using  opportunistic 
transmissions.  Being  a  passive  system  allows  the  PBR  receiver  to  remain  covert,  making 
it  more  resilient  to  detection  and  EA.  The  ability  to  leverage  available  transmission  and  to 
detect  targets  passively  also  serves  to  enhance  situational  awareness,  thereby  advancing 
one’s  position  in  the  OODA  loop  during  an  operation.  It  is  also  advantageous  to  use 
multiple  transmitters  at  separate  locations  for  detection  and  target  location  estimation  as  it 
adds  spatial  diversity  and  reduces  the  effect  of  DOP. 

1.  Findings  for  Detection  Coverage 

Target  detection  performance  varies  with  the  number  of  transmitters,  their 
parameters,  and  receiver  parameters.  It  also  depends  upon  the  bistatic  RCS  scattering 
characteristic  of  the  target.  Detection  performance — received  target  echo  SNR  and  area 
with  good  target  echo  SNR — improves  with  the  number  of  opportunistic  transmitters  and 
is  observed  to  vary  significantly  with  a  slight  change  in  receive  angle  due  to  the  rapid 
variation  in  RCS  returns  versus  angle  from  a  target.  Overall  detection  coverage  area 
increases  as  transmitter(s)  to  receiver  range  increase  at  the  expense  of  overall  sensitivity. 
Detection  coverage  is  generally  the  best  when  the  target  is  located  within  the 
transmitters-receiver  cluster  with  the  highest  SNR  observed  when  the  transmitter-target- 
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receiver  is  in  the  forward-scattering  configuration.  Detection  gaps — defined  as  areas  with 
the  low  SNR  at  the  receiver — occur  when  the  target  is  outside  the  transmitters-receiver 
cluster  or  away  from  bearings  extending  from  the  transmitter-receiver  baseline.  Outside 
the  cluster  of  transmitters-receiver,  detection  is  better  when  the  transmitter-target-receiver 
is  in  the  back-scattering  configuration.  Other  than  the  areas  that  satisfy  a  back-scattering 
configuration,  it  is  observed  that  there  are  additional  areas  or  hyperbolic-like  curves — 
typical  of  TDOA  measurements — extending  from  the  transmitters  that  exhibit  good 
detection  coverage.  Therefore,  the  optimal  transmitters-target-receiver  configuration  and 
geometry  for  target  detection  coverage  is  such  that  the  target  is  within  the  spatial  “net”  of 
transmitters-receiver.  The  findings  with  regard  to  detection  coverage  align  with  the 
findings  and  observations  in  [4,  40]. 

2.  Findings  for  Target  Location 

Two  target  location  estimation  methods  were  investigated:  (1)  elliptical  method 
and  (2)  hyperbolic  method.  The  elliptical  method  solves  for  the  bistatic  triangle 
parameters  for  each  PBR  pair;  these  parameters  are  then  used  to  define  transmitter-to- 
target  and  receiver-to-target  bearings,  after  which  the  “best  fit”  point  of  intersection  is 
estimated  as  the  target  position.  In  the  hyperbolic  method,  time  delay  measurements  from 
all  PBR  pairs  are  used  to  estimate  the  target’s  location  by  finding  the  least-squares 
solution.  The  results  from  the  target  location  model  showed  that  the  DF  measurement 
error  is  the  dominant  source  of  error  for  the  elliptical  method,  while  the  time  delay 
measurement  error  is  the  primary  source  of  error  for  the  hyperbolic  method.  The 
hyperbolic  method  generally  estimates  the  target’s  location  more  accurately  than  the 
elliptical  method  as  it  uses  fewer  measurement  sources.  It  was  also  observed  that  overall 
location  accuracy  is  poor  when  measurement  errors  are  modeled  as  dependent  on  the 
SNR  of  the  received  target  echo. 

When  SNR-dependent  measurement  errors  are  considered,  results  from  the  target 
location  model  agree  closely  with  the  observations  made  from  the  detection  coverage 
results.  Target  location  methods  produce  the  most  accurate  results  when  the  target  lies 
within  the  transmitters-receiver  cluster  in  a  receiver-centered  geometry.  The  transmitters- 
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receiver  cluster  area  can  be  expanded  by  increasing  the  transmitters  to  receiver  range; 
however,  this  reduces  the  overall  location  accuracy  within  the  cluster. 

Monte  Carlo  simulations  are  performed  for  both  target  location  methods  to 
explore  the  relation  between  spatial  error  distribution  and  transmitter-target-receiver 
geometry.  The  results  from  Monte  Carlo  simulations  in  a  scatter  plot  reveal  the 
uncertainty  area  associated  with  the  target  location  estimate  for  a  particular  geometry. 
Across  all  geometries,  the  uncertainty  area  associated  with  estimates  from  the  hyperbolic 
method  covers  a  smaller  area  compared  with  that  from  the  elliptical  method,  suggesting 
that  the  hyperbolic  method  estimates  the  target  location  with  higher  accuracy.  For  targets 
within  the  cluster  of  transmitters-receiver,  the  uncertainty  area  would  be  approximately 
circular.  However,  when  targets  are  outside  the  cluster  of  transmitters-receiver,  the 
uncertainty  area  is  elliptical.  Graphical  illustrations  were  presented  to  understand  how 
measurement  errors  translate  to  uncertainty  estimates  for  both  methods.  Monte  Carlo 
simulation  results  from  the  target  location  model  showed  that  the  elliptical  and  hyperbolic 
target  location  methods  have  approximately  orthogonal  error  distribution,  which  are  more 
pronounced  when  the  uncertainty  area  associated  with  the  estimates  are  elliptical.  That  is, 
the  semi -major  axis  of  the  error  ellipse  of  estimates  from  the  elliptical  method  is 
approximately  aligned  with  the  semi-minor  axis  of  the  error  ellipse  of  estimates  from  the 
hyperbolic  method,  and  vice  versa. 

The  method  to  compute  the  error  ellipse  parameters  associated  with  the 
hyperbolic  method  was  presented  and  included  in  the  target  location  model.  Every  target 
position  estimated  from  the  hyperbolic  method  would  be  accompanied  by  a  measure  of 
confidence  in  the  form  of  an  error  ellipse.  Like  Monte  Carlo  simulations,  the  error  ellipse 
shows  how  measurement  errors  translate  to  an  estimate’s  uncertainty.  For  an  error  ellipse 
at  (l-a)xlOO%  confidence,  (1  -a)  of  the  100  Monte  Carlo  simulated  estimates  will,  on 
average,  fall  within  the  computed  ellipse. 

B.  FUTURE  WORK 

Detection  coverage  models  and  target  location  models  developed  in  this  thesis  can 
be  further  improved  to  include  bistatic  beam-to-beam  issues  and  complex  sources  of 
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measurement  and  data  error  such  as  sea  clutter  and  errors  associated  with  sea-borne 
vessel’s  motion  (e.g.,  pitch  and  roll).  Real  data  can  be  used  in  the  models  developed  in 
this  study  to  extend  the  findings  regarding  target  detection  coverage  and  location 
accuracy.  An  extension  to  the  current  work  is  to  investigate  the  feasibility  and 
performance  of  using  multiple  passive  cooperative  opportunistic  transmitters  in  maritime 
target  detection  and  location  as  proposed  in  [40],  A  combined  target  location  method  that 
leverages  on  the  spatial  precision  of  the  elliptical  and  hyperbolic  target  location  methods 
can  be  explored  to  improve  estimation  accuracy.  Alternatively,  the  use  of  other  sources  of 
measurements,  like  Doppler  information  to  estimate  the  target’s  location,  can  be 
investigated. 

With  target  detection  and  location  in  place,  the  next  step  is  to  explore  tracking 
algorithms,  classification  methods,  and  identification  methods.  The  MIMO  tracking 
algorithm,  which  uses  the  coherent-on-receiver  method  [41]  or  the  track-before-detect 
method  that  enables  tracking  and  detection  of  targets  in  low  SNR  environments  [42],  are 
feasible  options.  As  stated  in  [4],  an  interesting  extension  would  be  to  consider  the  use  of 
bistatic  ISAR  signatures  in  the  forward-scattering  and  back-scattering  configurations  to 
aid  in  target  identification  and  classification  [25], 
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APPENDIX  A.  DERIVATION  OF  ERROR  ELLIPSE  PARAMETERS 
FROM  BIVARIATE  NORMAL  DISTRIBUTION 


This  appendix  presents  the  derivation  of  parameters  that  define  the  error  ellipse 
[38].  To  understand  the  how  measurement  errors  translate  to  estimate  position  uncertainty 
in  the  form  of  an  error  ellipse,  consider  a  univariate  normal  distribution  with  mean  p  and 

variance  cr2  for  the  random  variable  X  : 


f(x)  = 


1 


x/2  ~na2 


(X-JLI) 

2<j2 


(66) 


for  -oo  <  x  <  go  ,  -oo  <  ju<  oo  and  cr2  >  0 . 

Similarly,  a  multivariate  normal  distribution  with  mean  p  and  covariance  S  for 
random  variables  in  X  is  expressed  as 


/XX)  = 


1 

- e  2 


yl@n)p  p 


(67) 


for  — OO  <  X  <  oo ,  —oo  <  p  <  x  and  for  positive  definite  Z  •  Eq.  67  is  derived  from  Eq.  66 
by  rewriting  the  term  ixX]  in  linear  algebra  formulation: 

(x-juf  (c>2)  \x-ju). 


(x-p)rz-‘(x-p). 


(68) 

(69) 


Eq.  69  is  the  squared  Mahalanobis  distance  where  X  and  p  are  (pxl)  vectors, 


Z  is  a  (p~xp)  positive  definite  covariance  matrix,  and  p  denotes  the  number  of 
variables  in  vector  X . 


To  derive  the  semi -major  axis  and  semi -minor  axis  lengths  of  an  error  ellipse  in 
the  XY-plane,  reduce  Eq.  67  into  a  bivariate  normal  distribution  where 


X 

,1  = 

Pxy^x^y 

_y_ 

_My_ 

paa 

•  yx  y  x 

J 

X  = 


By  substituting  Eq.  70  into  the  Mahalanobis  distance  in  Eq.  69, 


(70) 
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(x-n)  s-,(x-n)=  — 


p 


xy 


(x-pxf  2pxy(x-jux)(x-juy)  (. x-jUy ) 


+  - 


<T  CT 

X  J 


(71) 


the  probability  density  function  of  a  pair  of  jointly  Gaussian  random  variables  is  then  expressed 
as 


fxy(x>y) 


(_  I 

(. X~JUX  f  2  (  X~JUx  )(X-JUy  )  1  (  X-jUy  ^ 

[  2(1  -p%) 

Of  ^x^.)  O y 

for  —oo  <  x  <  oo  ,  — oo  <  y  <  oo . 


(72) 


Figure  49  shows  the  joint  Gaussian  probability  density  function  (pdf)  for  various 
correlation  and  standard  deviation  values.  The  shape  and  orientation  of  the  joint  Gaussian 
pdf  surface  depends  on  the  values  of  o x ,  <y  r ,  and  pxy ,  while  the  size  of  the  surface 

contours  depends  on  the  pdf  value  fXY  (x,  y)  in  consideration. 


An  error  ellipse  is  defined  as  a  locus  of  points  with  constant  /’yy(x,>’)  values,  as 
seen  in  the  surface  contour  plots  in  Figure  49.  Graphically,  the  error  ellipse  is  the  points 
of  intersection  between  the  pdf  surface  and  a  horizontal  plane.  Given  the  values  of  ox , 

<jy ,  and  p ^ ,  (xy  (-A  y)  is  constant  for  values  of  x  and  y  when  the  exponential  argument 
in  Eq.  72  is  constant.  Mathematically,  this  is  expressed  as 

(x-^)rz-'(x-t.)  =  c.  (73) 

Let  Xj  =  X-fi  such  that  Eq.  Error!  Reference  source  not  found,  reduces  to 

X[  z_1  x,  =  c. 

Factorize  the  covariance  matrix  X  by  eigen-decomposition: 


(74) 


X  =  PDP  , 


(75) 


where  P  is  a  matrix  of  eigenvectors  and  D  is  a  diagonal  matrix  whose  elements  are  the 
corresponding  eigenvalues  (/lj ,  ) . 
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Figure  49.  Joint  Gaussian  pdf  surface  and  contours  for  various  <jx  ,  oy ,  and  pxv 

values.  After  [38]. 
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The  eigenvectors  of  the  covariance  matrix  represent  the  direction  in  which  the 
data  varies  the  most  and  are  the  vectors  along  the  error  ellipse’s  semi-major  and  semi¬ 
minor  axis.  The  corresponding  eigenvalues  indicate  the  spread  of  the  data  in  the  direction 
of  the  eigenvectors  (Figure  50).  For  an  axis-aligned  error  ellipse  such  that  covariance  is 
zero,  the  eigenvalues  equal  the  variances  of  the  covariance  matrix,  and  the  eigenvectors 
are  in  the  direction  of  the  x-axis  and  y-axis.  However,  for  correlated  data,  the 
eigenvectors  represent  the  direction  with  the  largest  spread  of  data  and  the  eigenvalue 
determines  the  extend  of  spread. 


Figure  50.  Eigenvectors  (A1,A7)  of  a  covariance  matrix  on  an  error  ellipse. 


If  X  is  a  symmetric  matrix,  P  1  =  Pr,  Eq.  75  is  then  expressed  as 

Z  =  PDPr.  (76) 

When  X  can  be  eigen-decomposed  and  when  none  of  its  eigenvalues  are  zero,  then  X  is 
non-singular  and  its  inverse  is  given  by 

S-'=PD-'Pr  (77) 

and 

DWI-P.  (78) 
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As  illustrated  in  Figure  50,  define  the  r-axis  and  5-axis  as  the  principle  axis  of  the 
error  ellipse.  The  principle  axis  is  defined  as  rotating  the  original  x-axis  and  y-axis 
through  an  angle  6  and  is  expressed  mathematically  as 


=  P 


(79) 

To  determine  the  semi-major  and  semi -minor  axis  lengths  of  the  error  ellipse, 
substitute  Eq.  78  and  Eq.  79  into  Eq.  74 

[r  5]Pr  Z_1  P[/ 

[r  s]D~l[r 


\r  s 


r- 


s]T  =c 


[r  s] 


% 


a 


[r  sf  =c 


r  s 
- 1 - =  c. 


(80) 


The  general  form  of  an  ellipse  centered  at  the  origin  is 


2  2 

X  7  , 

- f  —  =  1 

a1  b2  (81) 

where  a  and  b  are  the  semi-major  and  semi-minor  axis  lengths,  respectively.  Comparing 
Eq.  80  and  Eq.  81, 

a  =  ^[c\\ 

b  =  Jcl  I 

V  Ai>A2  (82) 

To  understand  the  significance  of  the  constant  c ,  consider  a  (px  1)  multivariate 
normal  vector  X  ~  iV(ji,  Z) .  Geometrically,  degrees  of  freedom  is  interpreted  as  the 
dimension  of  the  vector  subspace.  Hence,  the  squared  Mahalanobis  distance  between  X 
and  n  is  a  chi-square  distribution  with  p  degrees  of  freedom  (Figure  51): 

(X-n)rI-'(X-n)~^.  (g3) 
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Chi-Square  PDF  with  p  Degrees  of  Freedom 


Figure  51.  Chi-square  pdf  for  p  degrees  of  freedom. 


For  a  hyper-ellipse,  defined  as  the  squared  Mahalanobis  distance  equal  to  a 
critical  value  of  a  chi-square  distribution  with  p  degrees  of  freedom  and  evaluated  at  a  , 
the  probability  of  the  random  vector  X  falling  inside  the  ellipse  is  1  —  a . 


Prj(X-(i)rI-|(X-(i)<^[  =  l- 


a 


(84) 


where  a  refers  to  the  area  to  the  right  of  the  chi-square  critical  value  (Figure  52). 

Therefore,  the  hyper-ellipse  is  defined  as  a  locus  of  points  with  constant 
fXY  (x,  y)  values  such  that 

(x-(l)rZ-,(x-(i)  =  4.  (85) 

By  comparing  Eq.  73  and  Eq.  85,  Eq.  82  can  be  written  as 

a  = 

> 

b  =  sjX 2,a^2  .  . 

y  (86) 

For  a  two-dimensional  hyper-ellipse,  the  chi-squared  distribution  with  two 
degrees  of  freedom,  p  =  2  will  be  used. 
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Chi-Square  PDF  with  2  Degrees  of  Freedom 


Figure  52.  Chi-square  pdf  where  the  area  to  the  right  of  the  critical  value  is  a  . 

For  an  error  ellipse  with  confidence  level  (1  -  a)  x  \  00% ,  the  relevant  x\M  value 

in  Eq.  86  is  obtained  from  the  chi-square  distribution  table.  A  sample  of  a  chi-square 
distribution  table  is  given  in  Table  8. 


Table  8.  Chi-square  distribution  table 


Degrees  of 
Freedom 

Probability 

0.95 

0.90 

0.80 

0.70 

0.50 

0.30 

0.20 

0.10 

0.05 

0.01 

0.001 

1 

0.004 

0.02 

0.06 

0.15 

0.46 

1.07 

1.64 

2.71 

3.84 

6.64 

10.83 

2 

0.10 

0.21 

0.45 

0.71 

1.39 

2.41 

3.22 

4.60 

5.99 

9.21 

13.82 

3 

0.35 

0.58 

1.01 

1.42 

2.37 

3.66 

4.64 

6.25 

7.82 

11.34 

16.27 

4 

0.71 

1.06 

1.65 

2.20 

3.36 

4.88 

5.99 

7.78 

9.49 

13.28 

18.47 

5 

1.14 

1.61 

2.34 

3.00 

4.35 

6.06 

7.29 

9.24 

11.07 

15.09 

20.52 

Referring  to  Figure  53,  the  error  ellipse’s  principle  axis  [r  s]  is  defined  as  the 

counterclockwise  rotation  of  the  original  xy-axis  [x  v]  through  an  angle  9 .  The  vector 

representation  of  the  original  .^-coordinates  as  a  function  of  the  /^'-coordinates  is 

q  =  Ay  (87) 
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where  q  =  [r 


s]r ,  y  =  [x  yf  and  A 


cos  #  sin  # 
-sin#  cos# 


Figure  53.  Rotation  of  error  ellipse  principle  axis.  From  [38], 


From  Eq.  87, 

=  AZy  Ar, 


O  POO 

Z  x  r  xy  x  y 

y  2 

poo  o 

_r  xy  x  y  y 

Substituting  Eq.  89  and  Eq.  90  into  Eq.  88  yields 


a 

o 

_ 1 

cos# 

sin# 

poo 

r  xy  x  y 

cos# 

-sin# 

1 - 

o 

•a, 

1 _ 

-sin# 

cos# 

Pxy^x^y 

J 

sin# 

cos# 

(88) 

(89) 

(90) 
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which  is  further  reduced  to 


0  0  0  0  0 

<jr  =  ax  cos  0  +  2 pxvaxav  sin  0cos0  +  ov  sin  0, 


0  0  0  0  0 

<rs  =  <jx  sin  0-2 p  <rxi 7y  sin 0cos0- cry  cos  0, 


0  =  (&y  ~  ^  )  sin  0 cos  0 + pxy(Jx(Jy  (cos2  0  -  sin2  0). 


Simplifying  Eq.  92  gives  the  orientation  of  the  error  ellipse: 


#  =  -tan  1 
2 


^  2p  <j  <7  ^ 

^  xeye  xe  y, 


V 


(91) 

(92) 

(93) 


(94) 
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APPENDIX  B.  SPECIFICATION  SHEET  FOR  MANTADIGITAL 
RADAR  BY  KELVIN  HUGHES 


NAVIGATION 

RADAR 


MANTADIGITAL™  RADAR  SYSTEMS 


van*r*  rrvO»  rf  25*^  r  «  t*nd  magn*«ror  rvl  £narpC/eTV  200VV  aort  ?tnin  vwr*o<**+ 

DISPLAY  TECHNOLOGY 


to}M  ttur  TtHi 
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C'VMrr>aci  Target  Cdtacxini  oa  aeaactcn  chan  w  ^anr  0*4 


ENHANCED  TARGET  DETECTION 

FT  r*  .  .  .  ^.rV-y  Ih. 

u:n*  *l<^j  nlj’tu  ia***  a  brf  |«toJ  *d  li'a  atd  P*iT  •lxaln|  iittur  il 
oiAUi*  »rf  il«*.n*r)  •  ■d  r-issm)  i mrjm%  AV«n  vwtTi  V»  <1#  PP  I 


•lardy  >y  «  mU»>  «^J  u* 
■-*"8  tt+rmx/to+mf 

^n3  yar^ng  nao^  t 


DUAL  PPI 

f  rtfi  r^odatn  radii  »scMv  "**»  daoonvng  cusirod  w*t\  riymscon  such  a  torn  ASP*  >f«woo»  t>  K  A1S  and  Oar  -mu  fl^iao  co*?rvj  Vv 

Id*  hull  oar  La  ail  lli  IVji  _/*  om  and  IViilm  f»  ofatr  ^wn)  Vw  a  tall#  cd  ndlffl  Aui  mnXilliy*  u/  Ma  trf 

a'»l  o <xr*a  and  dcmo  u»  am*  mu**. 

SHARPEYE™  SOLID  STATE  TECHNOLOGY 

Me  n^gr—  rror  rrnian  rq  jt««  rar»nrcw  rt'a  n*jn  and  on  coat  nf  wwardif. 

wmt^yyoy  ntsIMr  Th*  **gr*  «•  dtr**.  am  rac**  c*  »tKJ»ne  ninvirCT  hmw  *r»-  n?«  to  utwwdrttr  ac«rs  T**  manse****  •  »wi 


ANTENNA  ANO  TURNING  UNITS 


•adng  anr  rrnriM  *>  %arvca  da  Tne  j'^mw  i aa  ma  ro^aod  »ff~s  im 
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MANTADIGITAL™  RANGE 


Figure  54.  Technical  Specifications  for  MantaDigital  Radar  by  Kelvin  Hughes. 

From  [43], 
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APPENDIX  C.  SPECIFICATION  SHEET  FOR  TELEDYNE 
DEFENCE  QR026  EW  RECEIVER 
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TELEDYNE 

DEFENCE 

A  Teiedyne  Technologies  Company 


Figure  55.  Technical  Specifications  for  Teledyne  Defence  QR026  EW  Receiver. 

From  [44], 
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APPENDIX  D.  SPECIFICATION  SHEET  FOR  POYNTING 
DEFENCE  DF-A0062  DF  RECEIVER 


Poynting 


Antennas 


VERSION  2,0 


Direction  Finding  and 
Monitoring  Antenna 

20  -  6000  MHz _ 

Product  Code:  DF-A0062 


PRODUCT  DESCRIPTION: 
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Figure  56.  Technical  Specifications  for  Poynting  Defence  DF  A0062  DF 

Receiver.  From  [45], 
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APPENDIX  E.  DETECTION  COVERAGE  PLOTS  FOR  TARGET 

PLANE  IN  S-BAND 


Referring  to  detection  coverage  for  the  plate  at  0°  and  180°  orientations  in  Figure 
57  and  Figure  58,  respectively,  received  SNR  is  the  greatest  when  the  transmission  and 
reception  of  the  target  pane  echo  occurs  in  the  vicinity  of  the  target  plane  faces  (at  90° 
and  270°  RCS  mainlobes)  and  minimum  when  transmission  and  reception  occurs  at  the 
edges  of  the  plate  (0°  and  180°).  Similar  observations  were  made  for  detection  coverage 
for  target  plane  at  90°  and  -90°  orientation  (Figure  59  and  Figure  60). 

Comparing  the  detection  coverage  for  a  plate  at  a  45°  orientation  with  different 
transmitter  positions  in  Figure  61  and  Figure  62,  it  is  evident  that  different  transmitter- 
receiver  geometry  can  result  in  similar  detection  coverage.  The  target  position  with  the 
greatest  SNR  is  observed  when  either  transmission  or  reception  LOS  falls  in  the  vicinity 
of  the  plate  faces.  Similar  observations  were  made  for  a  plate  at  a  -45°  orientation  (Figure 
63  and  Figure  64). 


Figure  57.  Detection  coverage  (left)  for  Tx  at  [-5000,  0]  and  target  plate  at  0° 
orientation  with  corresponding  S-band  bistatic  RCS  (right). 
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figure  58. 


Detection  coverage  (left)  for  Tx  at  [-5000,  0]  and  target  plate  at  0° 
orientation  with  corresponding  S-band  bistatic  RCS  (right). 
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Figure  59.  Detection  coverage  (left)  for  Tx  at  [-5000,  0]  and  target  plate  at  90° 
orientation  with  corresponding  S-band  bistatic  RCS  (right). 
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Figure  60.  Detection  coverage  (left)  for  Tx  at  [-5000,  0]  and  target  plate  at  -90° 
orientation  with  corresponding  S-band  bistatic  RCS  (right). 
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Figure  61.  Detection  coverage  (left)  for  Tx  at  [-5000,  0]  and  target  plate  at  45° 
orientation  with  corresponding  S-band  bistatic  RCS  (right). 
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Figure  62.  Detection  coverage  (left)  for  Tx  at  [0,  5000]  and  target  plate  at  45° 
orientation  with  corresponding  S-band  bistatic  RCS  (right). 
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Figure  63.  Detection  coverage  (left)  for  Tx  at  [-5000,0]  and  target  plate  at  -45° 
orientation  with  corresponding  S-band  bistatic  RCS  (right). 
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Figure  64.  Detection  coverage  (left)  for  Tx  at  [0,  5000]  and  target  plane  at  -45° 
orientation  with  corresponding  S-band  bistatic  RCS  (right). 
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APPENDIX  F.  DETECTION  COVERAGE  RESULTS 


F.l  VARY  NUMBER  OF  TRANSMITTERS 
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Figure  65.  Detection  coverage  at  S-Band  for  1  Tx  at  5  km  range  and  target  at  0° 

orientation. 
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Figure  66.  Detection  coverage  at  S-Band  for  1  Tx  at  5  km  range  and  target  at 

-90°  orientation. 

99 


Detection  Coverage  SNR  «;dB> 


Figure  67.  Detection  coverage  at  X-Band  for  1  Tx  at  5  km  range  and  target  at  0° 

orientation. 
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Figure  68.  Detection  coverage  at  X-Band  for  1  Tx  at  5  km  range  and  target  at 

-90°  orientation. 


100 


Octane*  <m)  ,  ,9* 

Figure  69.  Detection  coverage  at  S-Band  for  2  Tx  at  5  km  range  and  target  at  0° 

orientation. 
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Figure  70.  Detection  coverage  at  S-Band  for  2  Tx  at  5  km  range  and  target  at 

-90°  orientation. 
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Figure  7 1 .  Detection  coverage  at  X-Band  for  2  Tx  at  5  km  range  and  target  at  0° 

orientation. 


Figure  72.  Detection  coverage  at  X-Band  for  2  Tx  at  5  km  range  and  target  at 

-90°  orientation. 


102 


CKlnnoe  <m)  ,  jq4 

Figure  73.  Detection  coverage  at  S-Band  for  3  Tx  at  5  km  range  and  target  at  0° 

orientation. 
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Figure  74.  Detection  coverage  at  S-Band  for  3  Tx  at  5  km  range  and  target  at 

-90°  orientation. 
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Figure  75.  Detection  coverage  at  X-Band  for  3  Tx  at  5  km  range  and  target  at  0° 

orientation. 
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Figure  76.  Detection  coverage  at  X-Band  for  3  Tx  at  5  km  range  and  target  at 

-90°  orientation. 
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Figure  77.  Detection  coverage  at  S-Band  for  4  Tx  at  5  km  range  and  target  at  0° 

orientation. 
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Figure  78.  Detection  coverage  at  S-Band  for  4  Tx  at  5  km  range  and  target  at 

-90°  orientation. 
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Figure  79.  Detection  coverage  at  X-Band  for  4  Tx  at  5  km  range  and  target  at  0° 

orientation. 
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Figure  80.  Detection  coverage  at  X-Band  for  4  Tx  at  5  km  range  and  target  at 

-90°  orientation. 
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Figure  8 1 .  Detection  coverage  at  S-Band  for  8  Tx  at  5  km  range  and  target  at  0° 

orientation. 
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Figure  82.  Detection  coverage  at  S-Band  for  8  Tx  at  5  km  range  and  target  at 

-90°  orientation. 
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Figure  83.  Detection  coverage  at  X-Band  for  8  Tx  at  5  km  range  and  target  at  0° 

orientation. 
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Figure  84.  Detection  coverage  at  X-Band  for  8  Tx  at  5  km  range  and  target  at 

-90°  orientation. 
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F.2  VARY  TRANSMITTER  RANGE 
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Figure  85.  Detection  coverage  at  S-Band  for  1  Tx  at  5  km  range  and  target  at  0° 

orientation. 
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Figure  86.  Detection  coverage  at  S-Band  for  1  Tx  at  15  km  range  and  target  at  0° 

orientation. 
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Figure  87.  Detection  coverage  at  S-Band  for  2  Tx  at  5  km  range  and  target  at  0° 

orientation. 
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Figure  88.  Detection  coverage  at  S-Band  for  2  Tx  at  15  km  range  and  target  at  0° 

orientation. 
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Figure  89.  Detection  coverage  at  S-Band  for  4  Tx  at  5  km  range  and  target  at  0° 

orientation. 


Figure  90.  Detection  coverage  at  S-Band  for  4  Tx  at  15  km  range  and  target  at 

0°  orientation. 
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Figure  91 .  Detection  coverage  at  S-Band  for  8  Tx  at  5  km  range  and  target  at  0° 

orientation. 


Figure  92.  Detection  coverage  at  S-Band  for  8  Tx  at  15  km  range  and  target  at  0° 

orientation. 
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F.3  VARY  TRANSMITTER-TARGET-RECEIVER  GEOMETRY 


Figure  93.  Detection  coverage  at  S-Band  for  4  Tx  in  receiver-centered  geometry 

and  target  at  0°  orientation. 


Figure  94.  Detection  coverage  at  S-Band  for  4  Tx  in  receiver-centered  geometry 

and  target  at  -90°  orientation. 
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Figure  95.  Detection  coverage  at  S-Band  for  4  Tx  in  receiver-centered  geometry 

and  target  at  45°  orientation. 


Figure  96.  Detection  coverage  at  S-Band  for  4  Tx  in  receiver-centered  geometry 

and  target  at  30°  orientation. 
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Figure  97.  Detection  coverage  at  S-Band  for  4  Tx  in  transmitter-clustered 

geometry  and  target  at  0°  orientation. 


Figure  98.  Detection  coverage  at  S-Band  for  4  Tx  in  transmitter-clustered 
geometry  and  target  at  -90°  orientation. 


115 


Figure  99.  Detection  coverage  at  S-Band  for  4  Tx  in  transmitter-clustered 

geometry  and  target  at  45°  orientation. 
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Figure  100.  Detection  coverage  at  S-Band  for  4  Tx  in  transmitter-clustered 

geometry  and  target  at  30°  orientation. 
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F.4  RANDOMLY  DISTRIBUTED  TRANSMITTERS 
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Figure  101.  Detection  coverage  at  S-Band  for  8  Tx  randomly  positioned  around 
Rx  at  5-20  km  range  and  target  at  0°  orientation. 
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Figure  102.  Detection  coverage  at  S-Band  for  8  Tx  randomly  positioned  around 
Rx  at  5-30  km  range  and  target  at  0°  orientation. 
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Figure  103.  Detection  coverage  at  S-Band  for  20  Tx  randomly  positioned  around 
Rx  at  5-20  km  range  and  target  at  0°  orientation. 
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Figure  104.  Detection  coverage  at  S-Band  for  20  Tx  randomly  positioned  around 
Rx  at  5-30  km  range  and  target  at  0°  orientation. 
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Figure  105.  Detection  coverage  at  S-Band  for  8  Tx  randomly  positioned  around 
Rx  at  5-20  km  range  and  target  at  -90°  orientation. 
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Figure  106.  Detection  coverage  at  S-Band  for  8  Tx  randomly  positioned  around 
Rx  at  5-30  km  range  and  target  at  -90°  orientation. 
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Figure  107.  Detection  coverage  at  S-Band  for  20  Tx  randomly  positioned  around 
Rx  at  5-20  km  range  and  target  at  -90°  orientation. 


Figure  108.  Detection  coverage  at  S-Band  for  20  Tx  randomly  positioned  around 
Rx  at  5-30  km  range  and  target  at  -90°  orientation. 
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Figure  109.  Detection  coverage  at  S-Band  for  8  Tx  randomly  positioned  around 
Rx  at  5-20  km  range  and  target  at  45°  orientation. 


Figure  110.  Detection  coverage  at  S-Band  for  8  Tx  randomly  positioned  around 
Rx  at  5-30  km  range  and  target  at  45°  orientation. 
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Figure  111.  Detection  coverage  at  S-Band  for  20  Tx  randomly  positioned  around 
Rx  at  5-20  km  range  and  target  at  45°  orientation. 
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Figure  1 12.  Detection  coverage  at  S-Band  for  20  Tx  randomly  positioned  around 
Rx  at  5-30  km  range  and  target  at  45°  orientation. 
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APPENDIX  G.  TARGET  PATH  DETECTION  RESULTS 
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Figure  1 13.  Detection  coverage  along  target’s  path  (2  km  resolution)  at  S-Band 

for  4  Tx.  RCS  computed  by  calling  FEKO. 
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Figure  1 14.  SNR  at  receiver  along  target’s  path  (2  km  resolution)  at  S-Band  for  4 

Tx.  RCS  computed  by  calling  FEKO. 
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Figure  115.  Detection  coverage  along  target’s  path  (100  m  resolution)  at  S-Band 
for  4  Tx.  RCS  extracted  from  pre-computed  RCS  table. 
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Figure  1 16.  SNR  at  receiver  along  target’s  path  (100  m  resolution)  at  S-Band  for 
4  Tx.  RCS  extracted  from  pre-computed  RCS  table. 
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Figure  117.  Detection  coverage  along  target’s  path  (100  m  resolution)  at  S-Band 
for  8  Tx  randomly  position.  RCS  extracted  from  pre-computed  RCS 

table. 
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Figure  118.  SNR  at  receiver  along  target’s  path  (100  m  resolution)  at  S-Band  for 
8  Tx  randomly  position.  RCS  extracted  from  pre-computed  RCS  table. 
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APPENDIX  H.  TARGET  LOCATION  RESULTS 
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Figure  1 19.  Target  location  errors  using  elliptical  method  (S-Band,  4  Tx  at  5  km 
range)  and  SNR-independent  measurement  errors. 
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Figure  120.  Target  location  errors  using  hyperbolic  method  (S-Band,  4  Tx  at  5 
km  range)  and  SNR-independent  measurement  errors. 
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Figure  121.  Target  location  errors  using  elliptical  method  (S-Band,  8  Tx  at  5  km 
range)  and  SNR-independent  measurement  errors. 
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Figure  122.  Target  location  errors  using  hyperbolic  method  (S-Band,  8  Tx  at  5 
km  range)  and  SNR-independent  measurement  errors. 
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Figure  123.  Target  location  errors  using  elliptical  method  (S-Band,  4  Tx  at  15 
km  range)  and  SNR-independent  measurement  errors. 
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Figure  124.  Target  location  errors  using  hyperbolic  method  (S-Band,  4  Tx  at  15 
km  range)  and  SNR-independent  measurement  errors. 
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Figure  125.  Target  location  errors  using  elliptical  method  (S-Band,  8  Tx  at  15 
km  range)  and  SNR-independent  measurement  errors. 
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Figure  126.  Target  location  errors  using  hyperbolic  method  (S-Band,  8  Tx  at  15 
km  range)  and  SNR-independent  measurement  errors. 
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Figure  127.  Target  location  errors  using  elliptical  method  (S-Band,  4  Tx 
clustered)  and  SNR-independent  measurement  errors. 
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Figure  128.  Target  location  errors  using  hyperbolic  method  (S-Band,  4  Tx 
clustered)  and  SNR-independent  measurement  errors. 
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Figure  129.  Target  location  errors  using  elliptical  method  (S-Band,  8  Tx 
clustered)  and  SNR-independent  measurement  errors. 
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Figure  130.  Target  location  errors  using  hyperbolic  method  (S-Band,  8  Tx 
clustered)  and  SNR-independent  measurement  errors. 
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Target  location  errors  using  elliptical  method  (S-Band,  4  Tx  at  5  km 
range)  and  SNR-dependent  measurement  errors. 
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Target  location  errors  using  hyperbolic  method  (S-Band,  4  Tx  at  5 
km  range)  and  SNR-dependent  measurement  errors. 
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Figure  133.  Target  location  errors  using  elliptical  method  (S-Band,  8  Tx  at  5  km 
range)  and  SNR-dependent  measurement  errors. 
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Figure  134.  Target  location  errors  using  hyperbolic  method  (S-Band,  8  Tx  at  5 
km  range)  and  SNR-dependent  measurement  errors. 
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Figure  135.  Target  location  errors  using  elliptical  method  (S-Band,  4  Tx  at  15 
km  range)  and  SNR-dependent  measurement  errors. 
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Figure  136.  Target  location  errors  using  hyperbolic  method  (S-Band,  4  Tx  at  15 
km  range)  and  SNR-dependent  measurement  errors. 


135 


Target  Position  Estimate  Error  <m) 


Figure  137.  Target  location  errors  using  elliptical  method  (S-Band,  8  Tx  at  15 
km  range)  and  SNR-dependent  measurement  errors. 
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Figure  138.  Target  location  errors  using  hyperbolic  method  (S-Band,  8  Tx  at  15 
km  range)  and  SNR-dependent  measurement  errors. 
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T arpet  Posxion  Estimate  Error  (m) 
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Figure  139.  Target  location  errors  using  elliptical  method  (S-Band,  8  Tx  at 
random  positions)  and  SNR-dependent  measurement  errors. 
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Figure  140.  Target  location  errors  using  hyperbolic  method  (S-Band,  8  Tx  at 
random  positions)  and  SNR-dependent  measurement  errors. 
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Figure  141.  Target  location  errors  using  elliptical  method  (S-Band,  4  Tx 
clustered)  and  SNR-dependent  measurement  errors. 
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Figure  142.  Target  location  errors  using  hyperbolic  method  (S-Band,  4  Tx 
clustered)  and  SNR-dependent  measurement  errors. 
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Figure  143.  Target  location  errors  using  elliptical  method  (S-Band,  8  Tx 
clustered)  and  SNR-dependent  measurement  errors. 
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Figure  144.  Target  location  errors  using  hyperbolic  method  (S-Band,  8  Tx 
clustered)  and  SNR-dependent  measurement  errors. 
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